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Abstract. A review is given of the very recent developments in the fields of real and virtual Compton 
scattering off the nucleon. Both real and virtual Compton scattering reactions are discussed at low outgoing 
photon energy where one accesses polarizabilities of the nucleon. The real Compton scattering at large 
momentum transfer is discussed which is asymptotically a tool to obtain information on the valence quark 
wave function of the nucleon. The rapid developments in deeply virtual Compton scattering and associated 
meson electroproduction reactions at high energy, high photon virtuality and small momentum transfer 
to the nucleon are discussed. A unified theoretical description of those processes has emerged over the 
last few years, which gives access to new, generalized parton distributions. The experimental status and 
perspectives in these fields are also discussed. 



PACS. 13.60.Fz Elastic and Compton scattering 
12.38.Bx Perturbative calculations 



1 Introduction 

In the study of hadron structure, one of the main ques- 
tions is how hadrons and nuclei are built from quarks and 
gluons and how one goes over from a description in terms 
of quark and gluon degrees of freedom to a description in 
terms of hadronic degrees of freedom. 

Nowadays, precision experiments at high energy have 
established Quantum Chromo Dynamics (QCD) as the 
gauge theory of strong interactions describing the dynam- 
ics between colored quarks and gluons. QCD exhibits the 
property that the interaction between the quarks becomes 
weak at very short distances, which is known as asymp- 
totic freedom, and which allows us to use perturbation 
theory to describe high energy strong interaction phenom- 
ena. On the other hand at low energy, the QCD coupling 
constant grows, and quarks and gluons are confined into 
colorless mesons and baryons, which are the particles as 
seen through experiments. In this phase of hadronic mat- 
ter, the underlying chiral symmetry of QCD, due to the 
nearly massless up, down and approximately also strange 
quarks, is spontaneously broken. To describe this regime 
directly from the QCD Lagrangian is a hard task which 
still defies a solution due to the strong coupling constant 
requiring non-perturbative methods. Probably the most 
promising and direct approach is the numerical solution 
of QCD through lattice calculations. For static hadronic 
properties, such as e.g. masses, much progress has already 
been made, but for more complicated hadron structure 
quantities, such as e.g. nucleon parton distributions, lat- 
tice calculations are still in their infancy. 

In absence of a full numerical solution of QCD, which 
would be able to describe the rich complexity of the hadro- 
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nic many body systems from their underlying dynamics, 
a complementary strategy to refine our understanding of 
hadron structure is to perform new types of precision ex- 
periments in kinematical regimes at low energy, which re- 
quire an inherent non-perturbative description. Besides, 
one may perform new types of experiments at high en- 
ergies, in those kinematical regimes where factorization 
theorems have been proven. Such experiments will allow 
us to use an accurate perturbative QCD description of the 
reaction dynamics as a tool to extract new types of non- 
perturbative hadron structure information. 

In this quest at the intersection of particle and nuclear 
physics, the experiments performed with electromagnetic 
probes play an important role. A new generation of pre- 
cision experiments has become possible with the advent 
of the new electron accelerators and in combination with 
high precision and large acceptance detectors. In particu- 
lar, high precision Compton scattering experiments have 
become a reality in recent years. In Compton scattering, a 
real or virtual photon interacts with the nucleon and a real 
photon is emitted in the process. As this is a purely elec- 
tromagnetic process, it yields small cross sections (com- 
pared to hadronic reactions), but on the other hand con- 
stitutes a clean probe of hadron structure. 

In this paper, a review will be given of very recent 
developments in the field of real and virtual Compton 
scattering off the nucleon. I shall discuss real and vir- 
tual Compton scattering at the same time and point out 
their complementarity and the differences in the extracted 
nucleon structure information. The emphasis is on those 
kinematical regimes where a fruitful interpretation is terms 
of nucleon structure observables has been shown to be pos- 
sible. Virtual Compton scattering (VCS) off the nucleon 
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has been reviewed before in Ref. [Gui9£], which is referred 
to for most technical details. For the VCS part, the em- 
phasis is on the progress over the past two years in the 
light of the first high precision VCS data in the threshold 
regime now available, and on the rapid development in 
the field of deeply virtual Compton scattering (DVCS) at 
large photon virtuality. 

In section M, the real Compton scattering (RCS) pro- 
cess below pion threshold is discussed. In this regime, the 
RCS process can be interpreted as the global response of 
the nucleon to an applied electromagnetic field, which al- 
lows us to access global nucleon polarizabilities. A disper- 
sion relation formalism is described, which was developed 
as a method to minimize the model uncertainty in the ex- 
traction of nucleon polarizabilites from both unpolarized 
and polarized RCS data at low energy. 

In section [| the virtual Compton scattering (VCS) 
reaction at low energy is discussed. It consists of a gener- 
alization of RCS in which both energy and momentum of 
the virtual photon can be varied independently, which al- 
lows us to extract response functions, parametrized by the 
so-called generalized polarizabilities (CP's) of the nucleon. 
A first dedicated VCS experiment was performed at the 
MAMI accelerator, and two combinations of those GP's 
have been measured. Their values are compared with nu- 
cleon structure model predictions. Further experimental 
programs are underway at the major electron accelerators 
to measure the VCS observables. It is outlined how results 
of such experiments can be interpreted in terms of the nu- 
cleon GP's, and in particular how polarization observables 
can separate the different GP's. 

Besides the low energy region, RCS at high energy and 
large momentum transfer is a tool to access information 
on the partonic structure of the nucleon. In section |[ a 
leading order perturbative QCD calculation of RCS is de- 
scribed, which was developed to extract the valence quark 
wave function of the nucleon. Such an approach is com- 
pared with competing mechanisms, and it is pointed out 
how the planned experiments can shed light on this field. 

Section |^ discusses the recent developments in deeply 
virtual Compton scattering (DVCS) and associated meson 
electroproduction reactions at high energy, high photon 
virtuality Q 2 and small momentum transfer to the nu- 
cleon. It is shown how a unified theoretical description of 
those processes has recently emerged and how it gives ac- 
cess to new parton distributions, the so-called skewed par- 
ton distributions, which are generalizations of the usual 
parton distributions as known from inclusive deep inelas- 
tic lepton scattering. Leading order perturbative QCD cal- 
culations of DVCS and different meson electroproduction 
reactions, using an ansatz for the skewed parton distribu- 
tions, are discussed in the kinematical regimes accessible 
at present or planned facilities. The corrections to those 
leading order QCD amplitudes and further open questions 
in this field are touched on briefly. 

In section ^, the calculation of the QED radiative cor- 
rections to the VCS process is discussed, which is indis- 
pensable to accurately extract the nucleon structure in- 
formation from VCS experiments. 



Finally, the conclusions and perspectives are given in 
section 0. 



2 Real Compton scattering (RCS) and 
nucleon polarizabilities 

2.1 Introduction 

In the study of nucleon structure, real Compton scattering 
(RCS) at low energy is a precision tool to access global in- 
formation on the nucleon ground state and its excitation 
spectrum. RCS off the nucleon is determined by 6 inde- 
pendent helicity amplitudes, which are functions of two 
variables, e.g. the Lorentz invariant variables v (related to 
the lab energy of the incident photon) and t (related to the 
momentum transfer to the target). In the limit v — > 0, cor- 
responding to wavelengths much larger than the nucleon 
size, the general structure of these amplitudes is governed 
by low energy theorems (LET) based on Lorentz invari- 
ance, gauge invariance and crossing symmetry. These the- 
orems predict that the leading terms of an expansion in 
v are determined by global properties of the nucleon, i.e. 
its charge, mass and anomalous magnetic moment. Fur- 
thermore, the internal structure shows up only at relative 
order v 2 and can be parametrized in terms of the polariz- 
abilities. In this way, there appear 6 polarizabilities for the 
nucleon, the electric and magnetic (scalar) polarizabilities 
a and (3 respectively, familiar from classical physics, and 
4 spin (vector) polarizabilities 71 to 74, originating from 
the spin 1/2 nature of the nucleon. These polarizabilities 
describe the response of the nucleon's charge, magneti- 
zation, and spin distributions to an external quasistatic 
electromagnetic field. As such the polarizabilities are fun- 
damental structure constants of the composite system. 

The differential cross section for RCS in the limit v — > 
is given by the (model independent) Thomson term, as 
a consequence of the LET. The electric and magnetic po- 
larizabilities then appear, in a low-energy expansion for 
RCS, as interference between the Thomson term and the 
subleading terms, i.e. as contribution of 0{v 2 ) in the dif- 
ferential cross section. In this way, a and (3 can in principle 
be separated by studying the RCS angular distributions. 
However, it has never been possible to isolate this term 
and thus to determine the polarizabilities in a model in- 
dependent way. The obvious reason is that, for sufficiently 
small energies, say v < 40 MeV, the structure effects are 
extremely small and hence the statistical errors for the 
polarizabilities large. Therefore, one has to go to larger 
energies, where the higher terms in the expansion become 
increasingly important and where also the spin polarizabil- 
ities come into play. To determine the nucleon polarizabil- 
ities from RCS observables, a reliable estimate of higher 
terms in the energy is therefore of utmost importance. 
To this end, actual experiments were usually analyzed in 
an uns ubtracted dispersion relation formalism at fixed t 
[Lvo97|. Using such an analysis, the proton scalar polariz- 



abilities were derived from Compton scattering data below 
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pion threshold, with the results [Mac95| : 

a = (12.1 ± 0.8 ± 0.5) x 10~ 4 fm 3 
[3 = (2.1 + 0.8 T 0.5) x lO^fm 3 . 



(1) 



Very recently, new RCS data on the pr oton below pion 
threshold have become available [ Olm00| . These data in- 
crease the available world data set substantially, and yield, 
in an unsubtracted DR formalism, the results : 



a = (11.89 ± 0.57) x 10~ 4 fm 3 
P = (1.17 ± 0.75) x 10~ 4 fm 3 . 



(2) 



The sum of the scalar polarizabilities, which appears 
in the forward spin averaged Compton amplitude, can be 
determined directly from the total photoabsorption cross 
section by Baldin's sum rule [ Bal60 |, which yields a rather 
precise value : 



■0 



(14.2 ± 0.5) x 10~ 4 fm 3 , 
(13.69 ± 0.14) x 10~ 4 fm 3 



(3) 
(4) 



with (|) from Ref. jDam70| and (|) from Ref. flBabgg ]. 

Similarly, the proton forward spin polarizability can be 
evaluated by an integral over the difference of the absorp- 
tion cross sections in states with helicity 3/2 and 1/2, 



7o = 7i - 72 - 274 



-1.34 x 10~ 4 fm 4 
= -0.80 x 10~ 4 fm 4 



(5) 
(6) 



with (|) from Ref. [|5an94| and (f| from Ref. ||Dre00a |. 
While the predictions of Eqs. (^^) rely on pion photopro- 
duction multipoles, the helicity cross sections have now 
been directly determined at MAMI by scattering photon s 
with circular polarizations on polarized protons [Are00|. 
The contribution to 70 for the proton within the measured 
integration range (200 MeV < v < 800 MeV) is [ |Are00| : 



7o 



800 MeV 
200 MeV 



(-1.68 ±0.07) x 10~ 4 fm 4 



(7) 



The contribution below 200 MeV can be estimated with 



the HDT pion photoproduction multipoles [Han98| to yield 
+1.0, and an estimate of the contribution above 800 MeV 



based on the SAID pion photoproduction multipoles | Arn9( 
yields -0.02, which results in a total value : 70 = -0.7 (here 
and in the following, all spin polarizabilities are given in 
units 10- 4 fm 4 ). 

Furthermore, unpolarized RCS data in the Zi(1232)- 
region were used to give - within a dispersion relation for- 
malism - a first prediction for the so-called backward spin 
polarizability of the proton, i.e. the particular combina- 
tion 7^ = 71 + 72 + 2 74 ente ring the Compton spin-flip 
amplitude at = 180° jTon! 



1* = 



1 2 g 

27.1 ± 2.2(stat + syst) ' (model) 



x 10" 4 fm 4 



These values for the polarizabilities can be compared 
with our present day theoretical understanding from chi- 
ral perturbation theory (ChPT). A calculation to 0(p 4 ) 



in heavy baryon ChPT (HBChPT), where the expansion 
parameter p is an external momentum or the quark mass, 
yields (here and in the following, a and (3 are given in units 
10- 4 fm 3 ) : a = 10.5 ±2.0 and f3 = 3.5 ±3.6, the errors be- 
ing due to 4 counter terms entering to th at ord er, which 
were estimated by resonance saturation [Ber92|. One of 
these counter terms describes the large paramagnetic con- 
tribution of the Zi(1232) resonance, which is partly can- 
celled by large diamagnetic contributions of pion-nucleon 
(N7r)-loops. In view of the importance of the A resonance, 
the calculation was also done by including the A as a dy- 
namical degree of freedom. This adds a further expansion 
parameter, the difference of the A and nucleon masses ("e 
expansion"). A calculation to 0(e 3 ) yielded a = 12.2 + 
+ 4.2 = 16.4 and /3 = 1.2 + 7.2 + 0.7 = 9.1, the 3 sepa- 
rate terms referring to contributions of N7r-loops (which is 
the Q(p 3 ) result), Z\-pole terms, and Z\7r-loops | Hem97a , 
|Hom98j . These 0(e 3 ) predictions are clearly at variance 
with the data, in particular a + (3 — 25.5 is nearly twice 
the rather precise value determined from Baldin's sum rule 
Eq. %. 

The spin polarizabilities have also been calcul ated in 
HBChPT. The 0(e 3 ) predictions for the proton are fHem98fl 
70 = 4.6 - 2.4 - 0.2 + = +2.0, and 7Tr = 4.6 + 2.4 - 
0.2 — 43.5 = —36.7, the 4 separate contributions refer- 
ring to N7T-loops (0(p 3 ) result), Z\-poles, ZW-loops, and 
the triangle anomaly, in that order. It is obvious that the 
anomaly or 7r°-pole gives by far the most important contri- 
bution to 7^, and that it would require surprisingly large 
higher order contributions to bring 7Tr close to the value 
of Eq. (||) . Recently, the N7r-loop contribution to the spin 
polarizabilities ha yc been evaluated in HBCh PT to 0(p 4 ) 
by several groups |ji00| , |Vij00| , |Gel00| . In Refs. [jji00| , [Vij00f , 
the result for the proton is 70 = +4.5 — 8.4, where the 
two contributions are the 0(p 3 ) and 0(p 4 ) N7r-loop con- 
tributions, in this order. Bas ed on the large 0(p 4 ) cor- 
rection term, the authors in [ JiOO, Vij00| call the conver- 
gence of_the chiral expansion into question. However in 
Ref. iGelOfll , different results were obtained for the 0(p 4 ) 
N7r-loop contributions to the 4 spin polarizabilities. It was 
argued that these differences are due to how one defines 
and extracts the 0(p 4 ) spin-dependent polarizabilities in 
chira l effect ive field theories. Following the procedure of 
Ref. |GelO0H , which removes first all one-particle reducible 
contributions from the spin-dependent Compton ampli- 
tude, the resulting values for 70 and 7^ of the proton are 
7o = +4.6 - 5.6 = -1.0, and 77r = +4.6 - 1.2 = +3.4 
(without the 7r°-pole), the separate contributions being 
again the 0(p 3 ) and 0(p 4 ) N7r-loop contributions respec- 
tively. For 70, a c onverg ence of HBChPT at order 0(p 4 ) 
was not expected ]Gel00 |, whereas the result for 7^ - when 
adding the 7r°-pole contribution - is not compat ible with 
the estimate of Eq. M) obtained by Ref. [ TonS 



In order to refine our present understanding of the nu- 
cleon polarizabilities, a better understanding of the con- 
vergence of the HBChPT expansion is absolutely neces- 
sary, and it is to be hoped that a calculation to 0(e 4 ) will 
clarify the status. On the other hand, it is also indispens- 
able to minimize any model dependence in the extraction 
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of the polarizabilities from the data. To this end, a fixed-i 
subtracted di spersion relation (DR) formalism was devel- 
oped in Ref. | Drc00a | for RCS off the nucleon at photon 
energies below 500 MeV, as a formalism to extract the 
nucleon polarizabilities with a minimum of model depen- 
dence as is described in the following. 



2.2 Fixed-t subtracted dispersion relations for RCS 

To perform a dispersion theoretical analysis of Compton 
scattering, one has to calculate the 6 independe nt stru c- 
ture functions Ai{v, t), i = 1, 6 (defined in Ref. [Lvo97|). 
They are functions of the usual Mandelstam variable t, 
and of v, defined in terms of the Mandelstam variables 
s and u as v — (s — u) / '(4mjv), with m^r the nucleon 
mass. The invariant amplitudes Ai are free of kinematical 
singularities and constraints, and because of the crossing 
symmetry they satisfy the relation Ai(y,t) = Aj(— v, t). 
Assuming further analyticity and an appropriate high- 
energy behavior, the amplitudes Ai fulfill unsubtracted 
DR at fixed t : 



2 f 

ReAi(u,t) = Af(v,t) + -VI 



+ °° , , v' lm s Ai{v',t) 
dv' „ ' 



(9) 

where Af are the Born (nucleon pole) contributions, and 
where Im s Ai are the discontinuities across the s-channel 
cuts of the Compton process, starting from the thresh- 
old for pion production, vthr- However, such unsubtracted 
DR require that at high energies (y — -» oo) the ampli- 
tudes lm s Ai(v, t) drop fast enough such that the integral 
of Eq. (|^) is convergent and the contribution from the 
semi-circle at infinity can be neglected. For real Comp- 
ton scattering, Regge theory predicts t he follo wing high- 
energy behavior for v — > oo and fixed t |Lvo97| : 



*(*) 



A 



3,5,6 



*(t)-2 



*(*)-3 



/.' " . Ai ~ v 

(10) 

where a(t) < 1 is the Regge trajectory. Due to the high- 
energy behavior of Eq. Jig), the unsubtracted dispersion 
integral of Eq. (|J) diverges for the amplitudes A\ and A2 . 
In order to obtain useful results for these two amplitudes, 
L'vov et al. [ Lvo97| proposed to close the contour of the 
integral in Eq. (g) by a semi-circle of finite radius v max in 
the complex plane (instead of the usually assumed infinite 
radius!), i.e. the real parts of A\ and A2 are calculated 
from the decomposition 

ReA^t) = Af{y,t) + A$*{v,t) + A?<y t t) , (11) 

with A\ nt the s-channel integral from pion threshold v t hr 
to a finite upper limit u max , and an 'asymptotic contri- 
bution' Af s representing the contribution along the finite 
semi-circle of radius v max in the complex plane. In the ac- 
tual calculations, the s-channel integral is typically evalu- 
ated up to a maximum photon energy of about 1.5 GeV, 
for which the imaginary part of the amplitudes can be 
expressed through unitarity by meson photoproduction 
amplitudes (mainly l7r and 2ir photoproduction) taken 



from experiment. All contributions from higher energies 
are then absorbed in the asymptotic terms A1 S , which are 
replaced by a finite number of energy independent poles 
in the t channel. In particular the asymptotic part of A\ 
is parametrized by the exchange of a sc alar pa rticle in the 
t channel, i.e. an effective "cr meson" [Lvo97|. In a simi- 
lar way, the asymptotic part of A2 is described by the ir° 
f-channel pole. This procedure is relatively safe for A2 be- 
cause of the dominance of the ir° pole or triangle anomaly, 
which is well established both experimentally and on gen- 
eral grounds as Wess-Zumino-Witten term. However, it 
introduces a considerable model-dependence in the case 
of A x . 



It was therefore the aim of Ref. DreOOa | to avoid the 
convergence problem of unsubtracted DR and the phe- 
nomenology necessary to determine the asymptotic con- 
tribution. To this end, it was proposed to consider DR's 
at fixed t that are once subtracted at v = 0, 



ReAi(v, t) = Af{v,t) + [M0,t)-A?(0,t)] 

lm s Ai(i/,t) 



v 2 V 



dv' 



v 1 {v' 2 -v 2 ) 



(12) 



These subtracted DR should converge for all 6 invariant 
amplitudes due to the two additional powers of v' in the 
denominator, and they are essentially saturated by the irN 
intermediate states. In other words, the lesser known con- 
tributions of two and more pions as well as higher continua 
are small and may be treated reliably by simple models. 

The price to pay for this alternative is the appearance 
of the subtraction functions Ai(v = 0,i), which have to 
be determined at some small (negative) value of t. This 
was achieved by s etting u p once-subtracted DR, this time 
in the variable t [ DreOOa : 



i4 4 (0,i)-^f(0 > *) = a j + al~ pole 

(13) 

where the six coefficients a, = Ai(0, 0) — Af(0, 0) are sim- 
ply related to the six polarizabilities a, (3, 71, 72, 73, 74 (see 



Ref. I DreOOa I for details), and where a\ pole represents, in 



the case of A2, the contribution of the ir° pole in the t- 
channel. 

To evaluate the dispersion integrals, the imaginary part 
of the Compton amplitude due to the s-channel cuts in 
Eq. ( [l2|) is determined, through the unitarity relation, 
from the scattering amplitudes of photoproduction on the 
nucleon. Due to the energy denominator l/v'(v' 2 — v 2 ) in 
the subtracted dispersion integrals, the most important 
contribution is from the 7rA^ intermediate states, while 
mechanisms involving more pions or heavier meso ns in the 
intermediate states are largely suppressed. In Ref. jDrcOOa ] 
the ttN contribution was then evaluated usin g the p ion 
photoproduction multipolc amplitudes of Ref. [ Han98| | at 
photon energies below 500 MeV, and at the high er ener 
gies using the SAID multipoles (SP98K solution) [ Arn96 
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as input. The multipion channels (in particular the 7r7riV 
channels) were approximated by the inelastic decay chan- 
nels of the ttN resonances. It was found, however, that in 
the subtracted DR formalism, the sensitivity to the mul- 
tipion channels is very small and that subtracted DR are 
essentially saturated at v w 0.4 GeV. 

The subtracted i-channel dispersion integral in Eq. ( |l3| ) 
from 4m 2 to +00 is essentially saturated by the imaginary 
part of the t-channel amplitude 77 — * NN due to 7T7t inter- 
mediate states. The dependence of the subtraction func- 
tions on momentum transfer t can be calculated by includ- 
ing the experimental information on the i-channel process 
throu gh 7T7T i ntermediate states as 77 — > 7T7t — > NN. In 
Ref. [Drc00a|, a unitarized amplitude for the 77 — ► 7T7t 
subprocess was constructed, and a good description of the 
available data was found. This information is then com- 
bined with the 7T7T — > NN amplitudes determined from 
dispersion theory by analytical continuation of ttN scat- 
tering. In this way, one avoids the uncertainties in Comp- 
ton scattering associated with the two-pion continuum 
in the t channel, usually modeled through the exchange 
of a somewhat fictitious a meson. The second integral 



in Eq. (|13| ) extends from —00 to —2 (m^ + 2Mm n ) w 
—0.56 GeV 2 . As we address Compton scattering for pho- 
ton energies below about 500 MeV, the value of t stays 
sufficiently small so that the denominator in the integral 
provides a rather large suppression, resulting in a small 
contribution. The contributio n along t he negative i-cut is 
estimated in the calculations | Dre00b by saturation with 
A intermediate states. Altogether the remaining uncer- 
tainties in the s- and t- channel subtracted integrals due 
to unknown high-energy contributions, are estimated to 
be less than 1%. As a consequence, this subtracted DR 
formalism provides a direct cross check between Compton 
scattering and one-pion photoproduction. 

Although all 6 subtraction constants a\ to a§ of Eq. ( fi"3[ 
could be used as fit parameters in the present formal- 
ism, the fit was restricted to the par ameters a\ and 02, 
or equivalently to a — (3 and 7^ in [Dre00a|. The sub- 
traction constants 04, 05 and ae were calculated through 
an unsubtracted sum rule (Eq. (§) for v = t = 0). The 
remaining subtraction constant 03, related to a + (3 by 
a + — — (a .-; + a6)/(27r), was fixed through Baldin's sum 
rule |Bal60| , using the value a + (3 = 13.69 x 10~ 4 fm 3 
| Bab98T 



2.3 Results for RCS observables 

Since the polarizabilities enter as subtraction constants, 
the subtracted dispersion relations can be used to extract 
the nucleon polarizabilities from RCS data with a mini- 
mum of model dependence. The present formalism can be 
applied up to photon energies of about 500 MeV. 

Below pion production threshold, RCS data were usu- 
ally analyzed to extract a and (3. However, it was shown 
that the sensitivity to 7^ is not at all negligible, espe- 
cially at the backward angles and the higher energies, so 
that bot h a — (3 and 7^ should be fitted simultaneously 
[DreOOaj. 



RCS above pion threshold can serve as a complement 
to determine the polarizabilities, in particular the spin 
polarizabilities, and can provide valuable cross checks be- 
tween Compton scattering and pion photoproduction, pro- 
vided one can minimize the model uncertainties in the dis- 
persion formalism. The three types of dispersion integrals 
of Eqs. ( [l2] ) and ( |l3| ) in the formalism outlined here are 
evaluated as described above. As a representative result 
obtained within the subtracted DR formalism, the RCS 
differential cross sections above pion threshold are shown 
in Fig. [I] at fixed a — (3 = 10 (here and in the follow- 
ing in units of 10 -4 fm 3 ), while 7^ is varied between —27 
(here and in the following i n units of 10~ 4 fm 4 ) and —37 
(for more details, see Refs. |Dre00a , Drc00b|j ). By compar- 
ing all available data above pion threshold, it was con- 
cluded Drc00a| that there is no consistency between the 
pion photoproduction data from MAMI (entering through 
the dispersion integrals) and available Compton scatter- 
ing data, in particular when comparing with the LEGS 
data, which were used in the extraction of Eq. (^|) for 7^ . 
Therefore, new data in the A region a re called for, some of 
which have recently become available [ Wis99| . An analysis 
of those unpolarized data in a dispersion formalism favors 
also a much more negative value fo r th an extracted in 
Eq. (|). The fit performed by Ref. |Wis99| yields a-/3 = 
9.1 ± 1.7(stat + syst) ±1.2(mod), when using a value of 
7?r = -37.6. 



Y + p— >Y + p (a - p = 

1 70 1 



10) 




Fig. 1. Differential cross section for Compton scattering off 
the proton as function of the cm. photon angle for different 
lab energies. The total results of the subtracted DR formalism 
are presented for fixed a — f3 = 10 and different values of j w : 
j n — -37 (dashed-dotted lines), 7^ = -32 (full lines) and j w — 
-27 (das hed lines). The references to the data can be found in 
Dre00a |, 
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Besides the existing information from unpolarized data, 
a full study of the spin (or vector) polarizabilities will how- 
ever re quire dou ble polarization experiments. It was in fact 
shown [ DreOOa | that the scattering of polarized photons on 
polarized protons is very sensitive to 7^, in particular in 
the backward hemisphere and at energies between thresh- 
old and the A region. In addition, possible normalization 
problems can be avoided by measuring appropriate asym- 
metries. Therefore, future polarization experiments hold 
the promise to disentangle the scalar and vector polar- 
izabilities of the nucleon with the help of the described 
subtracted DR formalism, and to further quantify the nu- 
cleon spin response in an external electromagnetic field. 



2.4 Higher order polarizabilities of the proton 

As outlined above, the electric and magnetic polarizabili- 
ties arise as 0{v 2 ) corrections to the lowest order (Thom- 
son) scattering amplitude. One can then ask the question 
whether Compton scattering can also provide additional 
proton structure information via the use of higher-order 
polarizabilities. If one extends the a nalysis to consider 
spin-averaged 0(v A ) terms (see Ref. [HolOO for details), 
four new structures fulfill the requirements of gauge, P, 
and T invariance. Two new quantities, a p Ev and p Mlll rep- 
resent dispersive corrections to the lowest order static po- 
larizabilities, a and (3 respectively, and describe the re- 
sponse of the system to time-dependent fields. Two more 
quantities a p E2 , /3 P 42 , represent quadrupole polarizabilities 
and measure the electric and magnetic quadrupole mo- 
ments induced in a system by the presence of an applied 
field gradient. 

As to the experimental evaluation of such structure 
probes, it is, of course, in principle possible to extract 
them directly from Compton cross section measurements. 
However, it is already clear from the previous discussion of 
present data, that isolating such pieces from other terms 
which affect the cross section at energies above ~ 100 MeV 
is virtually impossible since additional higher order effects 
soon become equally important. Thus an alternative pro- 
cedure is required, which is made possible by the validity 
of dispersion relation s. Within the subtracted DR formal- 
ism of Ref. I DreOOa] outlined above, those higher order 
terms in the expansion of the Co mpton amplitudes Aj 
can be reasonably evaluated as in [|Hol00|| . These higher 
order polarizabilities can be expressed in terms of appro- 
priate derivatives of the RCS invariant amplitudes Aj at 
t,v 2 = 0, denoted by a^f, a,^ v . In terms of subtracted 
DR's, they take the form : 



ait 




dv 



,lm s A l (v',t = Q) 



(14) 



v 

-4Mm,-2 



dt' 



t' 2 



.(15) 



The subtracted DR in Eqs. (Jl4|, |l5|) were evaluated as de- 
scribed above and yield (all in units of 10~ 4 fm 5 ) : 

c? Ev = -3.84-0.19 + 0.06, 
0" Mv = +9.29 + 0.15 - 0.07, 



p 

X E2 — 



+29.31-0.10-0.17, 



fiP M2 = -24.33 + 0.10 - 0.34, 



(16) 



where the second and third entries on the rhs of Eq. ( |T^ ) 
estimate the uncertainties in the s- and i-channel disper- 
sion integrals. 



The values of Eq. (|16|) were then confronted in | HolOO 
with the predictions of HBChPT at 0(p 3 ) 



o( P 3 ) 



P 



l E2 



= 2.4, 
22.1, 



$^=3.7, 



'Mv 
Pi 



M2 



-9.5. 



(17) 



By comparing Eq. ( pL 6[ ) and (y_7() one finds that the size of 
a p E2 is about right, while for both (3 P M2 and (3 p Mll the sign 
and order of magnitude is correct but additional contri- 
butions are called for. The most serious problem lies in 
the experimental determination of which is negative 
in contradistinction to the chiral prediction and to sum 
rule arguments which assert its positivity. To see if inclu- 
sion of Zi(1232) degrees of freedom can help to resolve the 
ab ove dis crepancies, these quantities were also calculated 
in pJolOCH in HBChPT to 0(e 3 ), with the result : 



C(e 3 ) 



*-Eu 
p 
Y E2 



= 1.7, 
= 26.2, 



P Mv = 7-5, 
/& a = -12.3. 



(18) 



Except for the sign problem with a El/ indicated above, 
which persists in the e-expansion, the changes are gen- 
erally helpful, although the magnetic quadrupole polariz- 
ability is still somewhat underpredicted. 

In [ HolOO I, the described analysis was also extended to 
higher order contributions 0(v b ) to the proton spin po- 
larizabilities, for which 8 new structures were found. A 
dispersive evaluation of those higher order spin polariz- 
abilities showed a qualitative agreement with HBChPT 
(9(e 3 ) predictions. 

Recently, an evaluation of the higher order polarizabil- 
i ties of th e proton in HBChPT to 0(p 4 ) has been reported 
[HcmOOb], providing an important new testing g round for 
the chiral predictions. It was found [HcmOOt] that the 
(3(p 4 ) HBChPT result for the 4 quadrupole polarizabili- 
ties and the 8 spin polarizabilities at 0(v 5 ) of the proton 
are in e ncoura ging good agreement with the DR estimates 
of Ref. poTOpj . 

In summary, the subtracted DR formalism presented 
not only provides a formalism to extract the lowest or- 
der nucleon polarizabilities from present and forthcoming 
RCS data with a minimum of model dependence. It can 
also be used to obtain information about higher order po- 
larizabilities of the proton, in this way providing a great 
deal of additional nucleon structure information. 



The higher order polarizabilities are then obtained as lin - 
ear combinations of the a^t and a% tV (for details see [ Hol00|] ). 
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3 Virtual Compton scattering (VCS) and 
generalized nucleon polarizabilities 

3.1 Introduction 

The nucleon structure information obtained through RCS, 
as discussed in section |L can be generalized by virtual 
Compton scattering (VCS) below pion threshold. VCS can 
be interpreted as electron scattering off a target polarized 
by the presence of constant electric and magnetic fields. 
To see how VCS generalizes the RCS process, it is useful 
to think of the analogy with the electromagnetic form fac- 
tors. Their measurement through elastic electron-nucleon 
scattering reveals the spatial distribution of the charge 
and magnetization distributions of the target, whereas a 
real photon is only sensitive to the overall charge and mag- 
netization of the target. The physics addressed with VCS 
is then the same as if one were performing an elastic elec- 
tron scattering experiment on a target placed between the 
plates of a capacitor or between the poles of a magnet. In 
this way one studies the spatial distributions of the polar- 
ization densities of the target, by means of the generalized 
polarizabilities, which are functions of the square of the 
four- momentum, Q 2 , transferred by the electron. 

Experimentally, the VCS process is accessed through 
the electroproduction of photons, and we consider in all 
of the following the reaction on a proton target, i.e. the 
reaction ep — > epj. One immediately sees a difference 
with regard to the RCS jp — ► jp reaction, because in 
the ep — > ep7 reaction, the final photon can be emitted 
either by the proton, giving access to the VCS process, or 
by the electrons, which is referred to as the Bethe-Heitler 
(BH) process. The BH amplitude can be calculated ex- 
actly in QED, provided one knows the elastic form factors 
of the proton. Therefore it contains no new information on 
the structure. Unfortunately, light particles such as elec- 
trons radiate much more than the heavy proton. Therefore 
the BH process generally dominates or at least interferes 
strongly with the VCS process, and this may complicate 
the interpretation of the ep — ► epr/ reaction. The only way 
out of this problem is either to find kinematical regions 
where the BH process is suppressed or to have a very good 
theoretical control over the interference between the BH 
and the VCS amplitudes, as will be discussed below. 

Assuming that this problem is fixed, one can then pro- 
ceed to extract the nucleon structure information from 
VCS. In doing so, care has to be taken to separate the triv- 
ial response of the target, due to its global charge and/or 
a global magnetic moment. Indeed, if we put a proton 
in an electric field, the first effect we observe is that it 
moves as a whole. Similarly, the magnetic field produces 
a precession of the magnetic moment. This problem is ab- 
sent when one studies the polarizability of a macroscopic 
sample because it can be fixed in space by appropriate 
means, which is not possible for the proton. This absence 
of a restoring force explains why the trivial response due 
to the motion of charge and magnetic moment dominates 
over the response of the internal degrees of freedom. This 
is the p hysical origin of the low energy theorem (LET) 
| Low58| for VCS. All what is needed to calculate this part 



of the response, are the parameters which control the mo- 
tion, that is the mass, the charge, and the magnetic mo- 
ment. Once the motion is known, one can compute the 
amplitude for scattering an electron on this moving pro- 
ton, the so-called Born amplitude. Having separated the 
trivial response, one can parametrize the structure part of 
interest in the VCS process through t he so-ca lled general- 
ized polarizabilities (CP's) as in Ref. [Gui95|. 



3.2 Definition of generalized polarizabilities 

The known BH + Born parts of the VCS amplitude at low 
energy start at order 1 /q' in an expa nsion in the outgoing 
photon energy q'. The LET [Low58 asserts that the non 
trivial part of the VCS amplitude, the so-called non-Born 
part (denoted by Hnb), begins at order q'. There is of 
course also a contribution of order q' in the BH + Born 
amplitude, but this term is exactly known and therefore 
can be subtracted, at least in principle. So what is needed 
next is an adequate parametrization of -ff/vs ■ For this pur- 
pose, a mu ltipole expansion (in the cm. frame) was per- 
formed in [Gui95] in order to take advantage of angular 
momentum and parity conservation. The behaviour of the 
non-Born VCS amplitude Hnb at low energy (q' — > 0) 
but at arbitrary three-momentum q of the virtual photon, 
was then parametrized by 10 functions of q defined by : 



Limit of -J-H { /^ pL}S '-' 
q' q L NB 



(q', q) when q' 0. (19) 



In this notation, p (p 1 ) refers to the electric (2), mag- 
netic (1) or longitudinal (0) nature of the initial (final) 
photon, L (L 1 = 1) represents the angular momentum of 
the initial (final) photon, whereas S differentiates between 
the spin-flip (5 = 1) and non spin- flip (S — 0) charac- 
ter of the electromagnetic transition at the nucleon side. 
As the angular momentum of the outgoing photon is L' 
= 1, this leads to 10 q-dependent GP's, denoted generi- 
cally by P^ p L ' pL ^ s (q). By imposing the constraints due 
to nucleon crossing combined with charge conjugation in- 
variancc on the VCS amplitude, it was shown however in 
[ |Drc97t|Drc9"8b| th at 4 of t he GP's can be eliminated. Thus 
only 6 GP's , e.g. pui98| 



p(01,01)0 (q)) p(ii,U)o (q)) 

P (oi,oi)i (q)) P (iMi)i (q)j p(ib02)i (q)j p(01,12)l (q); (2Q) 

are necessary to give the low energy behaviour of Hjq b ■ 

In the limit q — > for the GP's, one finds the following 
relations with the polarizabilit ies (in gaussian units) of 
RCS, as discussed in section ^ [ Drc98c| : 



p(01,01)0/QN = 
p(H4DO (0) = 
p(01,12)l/ Q x = 
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P 



11,02)1 



(0) 



a, 



1 2^2 



(72 + 74) 



P (oi,oi)i (0)=0i 
P (n,H)i( 0)=0) 



(21) 



AMl xp (h,y)-AM* n+B °™(h,y) 
= 4(2h)K 2 [v v 1 ^/2e(l-e)Pt T (ci) + v\j\ - e 2 P^ T (q) 
■fvffVT^P^CcQ + sjV^H^fq)} , (24) 



where a em = 1/137.036 is the QED fine structure con- 
stant. 



3.3 VCS observables 

We next discuss how one can analyze ep — > epj observ- 
ables to extract the 6 GP's of Eq. @. 

The VCS unpolarized squared amplitude is denoted 
by M. Besides, one can consider VCS double polarization 
observables, which are denoted by AM(h,i) for an elec- 
tron of hclicity h, and which are defined as a difference 
of the squared amplitudes for recoil (or target) proton 
spin orientation in the direction and opposite to the axis 
i (i — x, y, z) (see Ref. | Vdh97a for details). In an expan- 
sion in q', M and AM. take the form 



M cxp 



M c * p 



AM cxp = 



q' 2 
AM C * P 



M c xp 



O(q'), 



AM e *l 
q' 



AM°* P + O(q') . (22) 



Due to the LET, the threshold coefficients M-2, M~\, 
AM -2, AM~i are known. The information on the GP's 
is contained in Mq XP and AMq* p . These coefficients con- 
tain a part which comes from the (BH+Born) amplitude 
and another one which is a linear combination of the GP's 
with coefficients determined by the kinematics. 

The unpolarized observable M^ xp was obtained by Ref. 
|Gui95 | in terms of 3 structure functions Pll{<1), Ptt(<1), 
Plt{<\), which are linear combinations of the 6 GP's, 



^BH+Born = 2R% [^(g) _ p TT ( q) ] 

3o 



+ Va - ™ V3 y/2e(l + e)P LT (q) , (23) 



M c xp 



where K2, e, qo, t>i, i>2, V3 are kinematical quantities (for 
details see Ref. pui98| ). 

The three double-polarizati on observ ables AM^ p (h, i) 
(i — x, y, z) were expressed by [ Vdh97a| in terms of three 
new independent structure functions P£ T (q), P L z T (q), and 
P i y(q), which are also linear combinations of the 6 GP's, 



AM cxp/u ~^ A *-< BH + Born 



l0 \h,z)-AM? 



(h,z) 



vT 



TT 



(q) + ^V / 2e(l-e)P£ T (q) 



+v 3 ^2e(l-e)P' L z T (q)} 



AM c xp (h,x) 



AM% H+BoTn (h,x) 



= A{2h)K 2 [vfy/2e(l-e)P£ T (q) + v% ^/l - £ 2 P^ T (q) 
+«f Vl~e 2 Ptt (q) + v% v/2e (1 - e) P'jjj. (q) } , 



where vf , v%, v 

other structure functions in Eq. (t 



1 ! •••! u 4 



are kinematical coefficients. The 



can be expressed 
^ pui98fl . There- 
fore, measuring those 6 structure functions amounts to 
determine the 6 independent GP's. 



in terms of Pll,Ptt, Plt, P[ t , P l z t , P lt 



3.4 Results for VCS observables below pion threshold 

In the previous sections, the observables of the ep — > ep'y 
reaction below pion threshold were outlined, and it was 
shown how the nucleon structure effect can be parametrized 
in terms of 6 independent GP's. 

To access the GP's, the experimental strategy of VCS 
in the threshold region consists of two steps. First, one 
measures the VCS cross section at several values of the 
outgoing photon energy. At low energies, the measurement 
of the VCS observables provides a test of the LET. Once 
the LET is verified, the relative effect of the GP's can be 
extracted using Eqs. ( p3| , p4[ ). 

The predictions for the Bethe-Heitler (BH) and Born 
cross sections below pion threshold are shown in Fig. 
The BH cross section has a characteristic angular shape 
and displays two "spikes" , which occur when the direction 
of the outgoing photon coincides with either the initial or 
final electron directions. In these regions, the cross section 
is completely dominated by the BH contributions. In or- 
der to determine the VCS contribution, one clearly has to 
minimize the BH contamination by detecting the photon 
in the half-plane opposite to the electron directions. 



P (e, e' p)y 




100 150 

cm, 

% (deg) 

Fig. 2. Results for the BH (dashed-dotted curve), Born 
(dashed curve) and BH+Born (full curve) p(e, e' p)7 differ- 
ential cross sections in MAMI kinematics : q = 600 MeV/c, 
q' = 111.5 MeV/c and e = 0.62, as function of the cm. an- 
gle 6>c.m. between real and virtual photon , and for in-plane 
kinematics. Calculations from Ref. [Vdh96]. 
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The first dedicated VCS experiment has been performed 
at MAMI, and for the first time two combinations (see 
Eq. (||)) of GP's have been determine d at Q 2 = 0.33 
GeV 2 and photon polarization e = 0.62 [RocOC], 



1 



Pll(Q ) - -Ptt{Q ) = (23.7 ± 2.2 ± 0.6 ± 4.3) GeV 



Plt(Q 2 



-5.0 ± 0.8 ± 1.1 ± 1.4) GeV 2 



(25) 



VCS ex periments at higher Q 2 (1-2 GeV 2 ) at JLab 



| Bert93fl , and at lower Q 2 at MIT-Bates J5ha97[ have al- 
ready been performed, and are under analysis at this time. 

The GP's have been calculated in various approaches 
and nucleon st ructure mode ls, ranging from constituent 
quark models [ Gui95 , Pas98 , a relativistic effective La- 



grangian model |Vdh96|, and the linear a- model Met96 
to ChPT |Hem97b| , pc'm00a[ . The GP's teach us about 
the interplay between nucleon-core excitations and pion- 
cloud effects, which are described differently in the various 
models. We focus here on the calculation of the GP's in 
HBChPT to 0(p 3 ), as it takes account of N7r-loop contri- 
butions in a systematic way. The 0(p 3 ) calculation yields 
for the two measured combinations at Q 2 = 0.33 GeV 2 
and e = 0.62 of Eq. @ the values jHem00a| : 



o( P 3 



P, 



LL 



Ptt/c 
Plt 



26.3 GeV , 
-5.7 GeV 2 



(26) 



which are in astonishing agreement with the experimen- 
tally determined values of Eq. (|2^). In particular, the 
0(p 3 ) ChPT calculation predicts quite large values for 
the spin GP's. As for the case of the RCS polarizabili- 
ties, the importance of the 0(p 4 ) corrections remains to 
be checked. 

If one wants to extract the different polarizabilitics 
from experiment, and in particular in case of the spin po- 
larizabilities, an unpolarized experiment is not sufficient 
as it gives access to 3 independent response functions only. 
To further separate the polarizabilities, one has to resort 
to double-polarization observables. Experimentally, at the 
existing high-duty-cycle electron facilities with polarized 
electron beams such as MAMI, MIT-Bates and JLab, dou- 
ble polarization VCS experiments can be performed by 
measuring the recoil polarization of the outgoing nucleon 
with a focal plane polarime ter. An e xperiment at MAMI 
has already been proposed [ d'Ho99b ], 

In Fig. H, the double-polarization asymmetry of Eq. ( |2^ ) 
with proton recoil polarization either along the x- or z- 
directions, is shown at q = 600 MeV/c for in-plane kine- 
matics. It is seen that the asymmetry yields a very large 
value (between 0.6 and 0.7) if the final proton is polarized 
parallel to the virtual photon. 



3.5 Dispersion relation formalism for VCS 

At present, VCS experiments at low outgoing photon ener- 
gies are analyzed in terms of low-energy expansions (LEXs) 
of Eq. (|22]). In the LEX, the non-Born response of the sys- 
tem to the quasi-constant electromagnetic field of the low 



> 

c 0. 



q = 600 MeV/c, q' = 

1 1 1 1 1 1 1 » » » 1 1 1 ^ 
p fm pol along z 



20 MeV/c, £ = 0.61 




Fig. 3. VCS double-polarization asymmetries (polarized elec- 
tron, proton recoil polarization along either the x- or z- 
directions) in MAMI kinematics (Q 2 = 0.33 GeV 2 ) as function 
of the cm. angle between real and virtual photon. The BH + 
Born result is shown by the dashed curves. The other curves 
show pr edictions from different model contributions calculated 
in Ref. |Vdh97a], to which we refer for details. 



energetic photon is proportional to the GP's, as expressed 
in Eqs. ( p3| , p4| ). As the sensitivity of the VCS cross sec- 
tions to the GP's grows with the photon energy, it is ad- 
vantageous to go to higher photon energies, provided one 
can keep the theoretical uncertainties under control when 
crossing the pion threshold. The situation can be com- 
pared to RCS, for which one uses a dispersion relation 
formalism as discussed in section ||, to extract the polar- 
izabilitics at energies above pion threshold, with generally 
larger effects on the observables. 

In order to go to higher energies and to check the valid- 
ity of LEXs at these higher energies, a dispersion r elation 



analysis for VCS has been developed very recently [PasOO 



Vdh00b |, which will allow to extract the GP's from data 
over a larger energy range. The same formalism also pro- 
vides for the first time a dispersive evaluation of 4 GP's. 

To set up a dispersion formalism for the VCS process, 
one starts from the helicity amplitudes : 



e 2 e L 



^(q,X)e;(q',X')u(p',s')M^u(p, S ), 

(27) 

with e the electric charge, q (q') the four-vectors of the 
virtual (real) photon in the VCS process, and p (p') the 
four-momenta of the initial (final) nucleons respectively. 
The nucleon hclicitics are denoted by s,s' = ±1/2, and 
u, u are the nucleon spinors. The initial virtual photon 
has helicity A = 0, ±1 and polarization vector e^, whereas 
the final real photon has helicity A' = ±1 and polarization 
vector e u . The VCS process is characterized by 12 inde- 
pendent helicity amplitudes T\i s >. >,«, . 

The VCS tensor Ad^ in Eq. (g7|) is then expanded into 
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a basis of 12 independent gauge invariant tensors 



i=l 



F 4 (QV,i)p: 



(28) 



as introduced in Ref. [Drc97| (starting from the ampli- 
tudes of Ref. p?ar75| ). The amplitudes F l in Eq. (g8j) 
contain all nucleon structure information and are func- 
tions of 3 invariants for the VCS process : Q 2 = —q 2 , 
v — (s — u)/(4mAr) which is odd under s «-> u crossing, 
and t. The Mandelstam invariants s, t and u for VCS are 



2m 2 N 



2 , and u = (q — p') 2 , 
— Q 2 , and tun is the 



defined by s = (q + p) 2 , t = (q 
with the constraint s + 1 + u 
nucleon mass. 

Nucleon crossing combined with charge conjugation 
provides the following constraints on the amplitudes Fi 
[] at arbitrary virtuality Q 2 : 

F l (Q 2 ,-v,t)=F t (Q 2 ,u 1 t) (i = l,...,12). (29) 



With the choice of the tensor basis of Ref. | Dre97 |, the 
resulting non-Born amplitudes F^ B (i = 1,...,12) are free 
of all kinematical singularities and constraints. 

Assuming further analyticity and an appropriate high- 
energy behavior, the non- Born amplitudes F^ B (Q 2 ,v, t) 
fulfill unsubtracted dispersion relations (DR's) with re- 
spect to the variable v at fixed t and fixed virtuality Q 2 : 



ReF l NB (Q 2 ,v,t) = -V 

TT 



dv 



,z/'Im s Fi(Q 2 ,i/',t) 



(30) 

with lm s Fi the discontinuities across the s-channel cuts 
of the VCS process. Since pion production is the first in- 
elastic channel, v t hr = m n + (m 2 + t/2 + Q 2 /2)/(2m^), 
where m w denotes the pion mass. 



The unsubtracted DR's of Eq. (30) require that at suf- 
ficiently high energies (y — > oo at fixed t and fixed Q 2 ) the 
amplitudes lm s Fi(Q 2 ,v,t) (i = 1,...,12) drop fast enough 
such that the integrals are convergent and the contribu- 
tions from the semi-circle at infinity can be neglected. It 
turns out that for two amplitudes, F\ and F$, an unsub- 
tracted dispersion integral as in Eq. (|30| ) does not exist 



| PasOO , whereas the other 10 amplitudes can be evaluated 
through unsubtracted dispersion integrals. This situation 
is similar as for RCS, where 2 of the 6 invariant amplitudes 
cannot be eva luated by unsubtracted dispersion relations 
either |Lvo97|| . 

The unsubtracted DR formalism for VCS also allows 
to predict 4 of the 6 GP's. The appropriate limit in the 
definition of the GP's is q' —* at finite q (see Eq. (|l9|)), 
which corresponds in terms of VCS invariants to v — > 
and t — > — Q 2 at finite Q 2 . One can therefore express the 
GP's in terms of the VCS amplitudes Fi at the point 
v = 0, t = —Q 2 at finite Q 2 , denoted in the follow- 
ing as : F(Q 2 ) = Ff B (Q 2 ,u = 0,t = -Q 2 ). The re- 
lations between the GP's and the Fi(Q 2 ) can be found in 



1 In JPasOOt , 4 of the 12 amplitudes of |5re97] were redefined 
by dividing them through v, such that all of them are even 
functions of v. This simplifies the formalism since only one 
type of dispersion integrals needs to be considered then. 



Ref. [ Dre97 . From the high-energy behavior for the VCS 
invariant amplitudes, it follows that one can evaluate the 
Fi (for i ^ 1, 5) through the unsubtracted DR's 



HQ 2 ) 



+ °V WWV.^-Q 2 ). (31) 

7 /' 



V 



Unsubtracted DR's for the GP's will therefore hold for 
those combinations of GP's that do not depend upon the 
amplitudes F\ and F5 n. Among the 6 GP's, the fo llowing 



4 combinations of GP's were found in Ref. [PasOO 



p(01,01)0 _j_ lp(ll,ll)0 _ 2 



E + m N \ 1/2 



+ <^F 2 + (2F 6 + F 9 
.% 



V3 V E 

)-Fl2 



m N q 



1/2 



p (oi,oi)i = _L( E + m N \ ^ 



3V2 

{{F 5 +F 7 + 4F n )+4m N F 12 } 



(32) 



(33) 



p(01,12)l 



1 



v/2 



-P 



(11,11)1 



go 



1 / E + m N 
3 V E 



1/2 



m N q 



x{(F 5 + F 7 + A Fn) + 4 m N (2 F 6 + Fg) } , (34) 

p (oi.i2)i , V3 p (n,02)i _]_( E + m N \ 1/2 go 
2 6\ E ) q 2 

x {qo (F 5 + F 7 + 4 Fu) + 8 m 2 N (2 F 6 + Fg) } , (35) 



where E = -\/q 2 + m 2 N denotes the initial proton cm. en- 
ergy, and qo = mpf — E the virtual photon cm. energy in 
the limit q' = 0. Unfortunately, the 4 combinations of GP's 
of Eqs. (H^)-(H^) can at present not yet be compared with 
the dat a. In particular, the only unpolarized experiment 
[|Roc00| measured two structure functions which cannot 
be evaluated in an unsubtracted DR formalism, as they 
contain in addition to p(oi.0i)0 + 1 / 2 p(n,ii)o of Eq ^ @ ) 

which is proportional to a + (3 at Q 2 = 0, also the gener- 
alization of a — p. 

The 4 combinations of GP's on the Ihs of Eqs. (|32])~ 
( [35|) can then be evaluated by unsubtracted DR's, from 
the dispersion integrals of Eq. ( J3l| ) for the Fi(Q 2 ). To this 
end, the imaginary parts lm s Fi in Eq. ( |3l| ) have to be 
calculated by use of unitarity. For the VCS helicity ampli- 
tudes of Eq. ( ^7|) (denoted for short by 7/ <), the unitarity 
equation reads : 



2 Im, T 



fi 



x 



(27r) 4 S i (P x -P i )Tj (f T x 



(36) 



where the sum runs over all po ssible intermediate states 
X that can be formed. In Ref. [PasOO , the dispersion in- 



tegrals of Eq. (|3l]) were saturated by the dominant contri- 
bution of the ttN intermediate states. For the pion photo- 
and clectroproduction helicity amplitudes in the range 



i*5 can appear however in the combination F5 + 4Fn, in 
which the 7r°-pole drops out, and which has a hi gh-ener gy be- 
havior leading to a convergent integral (see Ref. [ Pas00| ). 



Marc Vanderhaeghen: Real and Virtual Compton Scattering off the Nucleon 



11 



Q 2 < .5 GeV 2 , the phenomenological analysis of MAID 
[Drc99| was used, which contains both resonant and non- 
resonant pion production mechanisms. 



,101,01)1 




0.1 



p(0U2,l_ 1/( ^ p („,„ } , + Vj/2 p 



0.1 0.2 0.3 0.4 0.5 
(01,12)1 , ^ n 0(11,02)1 




Fig. 4. Dispersion results [ PasOC ] for 4 of the generalized po- 
larizabilities of the proton (full curves ) , compared with results 



of 0(p 3 ) HBCh PT |Hem97b| , [Hem00a| ] (dashed curves) and the 



linear a- model [ Mct96| (dashed-dotted curves 



In Figjj, the results for the 4 combinations of GP's 
of Eqs. (p2|)-(^5|) are shown in the DR formalism, and 
compared to the results of the Q( p 3 ) heavy-baryon chiral 
perturbation theory (HBC hPT) jHcm97b| , pcm00^ | and 
the linear cr-model [ Met96[ . The irN contribution to the 
sum p(Qi.oi)° + l/2P( 1 1 - 11 ) gives only about 80% of the 
Baldin sum rule | Bab98 |, because of a non- negligible high- 
energy contribution (of heavier intermediate states) to the 
photoabsorption cross section entering the sum rule, which 
is not estimated here. On the other hand, for the 3 com- 
binations of spin polarizabilities of Eqs. (j33|)-(|35|), the dis- 
persive estimates with irN states are expected to provide 
a rather reliable guidance. By comparing the DR results 
with those of HBChPT at 0(p 3 ), one remarks a rather 
good agreement for p(«M2)i + V^/2P (11 ' 02)1 , whereas for 
the GP's p(°i.oi)i and P<°M2)i _ l/(V2g )P (U ' n)1 , the 
dispersive results drop much faster with Q 2 . This trend is 
also seen in the relativistic linear cr-model, which takes 
account of some higher orders in the chiral expansion. 
It remains to be checked how the 0(p 4 ) corrections in 
HBChPT change this comparison with the DR estimates. 

To complete the DR formalism for VCS, one further- 
needs to construct the VCS amplitudes Pi and P5, for 
which the unsubtracted dispersion integrals of Eq. ( ^o|) 
do not converge. One strategy is to pr oceed i n an analo- 
gous way as has been proposed in Ref. [Lvo97| in the case 
of RCS. The unsubtracted dispersion integrals for Pi and 
P5 are evaluated along the real v-axis in a finite range 



gral along a semi-circle of finite radius v max in the complex 
^-plane is described by the asymptotic contribution f" s , 
which is parametrized by t-channel poles (e.g. for Q 2 = 0, 
Pf s corresponds to cr-exchange, and Fg s to 7r°-exchange). 

A full study of VCS observables within such a dis- 
persion formalism, including a parametrization o f the tw o 
asymptotic contributions, is presently underway |Dre00c|. 
This will yield a formalism to extract the nucleon GP's 
over a larger range of energies from both unpolarized and 
polarized VCS data. 



4 Compton scattering at large momentum 
transfer and the nucleon distribution 
amplitude 

4.1 Introduction 

Besides the low energy region, as discussed in section |^, 
RCS will also provide access to information on the par- 
tonic structure of the nucleon at sufficiently large momen- 
tum transfer. 

This regime is defined by requiring that all three Man- 
delstam variables (s, t, u) be large with respect to a typ- 
ical hadronic scale, say 1 GeV. In this case there is a 
prejudice ( actual ly proven in the case of elastic electron 
scattering [ 5te97 |) that the amplitude factorizes in a soft 
non-perturbative part, the distribution amplitude, and a 
hard scattering kernel which is calculable from perturba- 
tive QCD (PQCD). Because of asymptotic freedom, the 
perturbative approach must be to some degree relevant 
to the hard scattering regime. However, since the binding 
of the quarks and gluons in the hadrons is a long dis- 
tance, non-perturbative effect, the description of the reac- 
tion requires a consistent analysis of both large and small 
scales. When the reaction is hard enough, the relative ve- 
locities of the participating particles are nearly lightlikc. 
Time dilatation increases the lifetime of the quantum con- 
figurations which build the hadron. As a result, the par- 
tonic content, as seen by the other particles, is frozen. 
Moreover, due to the apparent contraction of the hadron 
size, the time during which momentum can be exchanged 
is decreased. Therefore one expects a lack of coherence 
between the long-distance confining effects and the short 
distance reaction. This incoherence between the soft and 
hard physics is the origin of the factorization which is il- 
lustrated in Fig. pL 




< v < +u n 



(with v n 



1.5 GeV). The inte- 



Long distance Short distance Long distance 

Fig. 5. Factorization of the RCS scattering amplitude in the 
hard scattering regime. 
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4.2 Factorization and the nucleon distribution 
amplitude 

The calculation of the RCS amplitude at large momentum 
transfer, follows the Brodsky-Lepage formalism [Bro80|, 
which leads to the factorized expression : 

T(X',h' N ,X,h N ) = 

dxidyj <j>* N (yj)T H (\' ,h' N ,yj,X,h N , af<; s, t) cf> N (xi), 

(37) 

where (xi, y{) are the momentum fractions of the quarks 
in the initial and final nucleon respectively, Th is the hard 
scattering kernel and <Pn is the distribution amplitude 
(DA) . The evaluation of Eq. (l37j ) requires a four- fold con- 
volution integral since there are two constraint equations 
(x% + x 2 + x 3 = 1 and 2/1 + 2/2 + 2/3 = 1)- In Eq. ( p7| ) a suf- 
ficiently large momentum transfer is assumed in order to 
neglect the transverse momentum dependence of the par- 
tons in the hard scattering amplitude Th- In this limit, 
the integration over the transverse momenta kj_^ (where 
kj^ = 0) acts only on the valence wavefunction 



W v (xi,x 2 ,x 3 ;'k ±1 ,'k_ L 2,'k± 3 ) , 



(38) 



which is the amplitude of the three quark state in the Fock 
expansion of the proton: 

\P> = &v\ qqq > +& q q I qqq, qq > I qqq, g > 

+... (39) 

This valence wave function Wy integrated up to a scale [i 
(which separates the soft and hard parts of the wavefunc- 
tion) defines the DA which appears in Eq. (pj 



^> N (Xi,fl) = 



d 2 k ±i <P v (x i ;'k ± i) 



(40) 



For ji much larger than the average value of the trans- 
verse momentum i n the p roton, this function 0n depends 
only weakly on /i [ Bro80| and this dependence can be ne- 
glected. 

The interest of the formalism is that the distribution 
amplitude is universal, that is independent of the partic- 
ular reaction considered. Several distributi on amplitudes 
have b een modeled using QCD sum rules HChc84phe89j 
Kin87]. They have a characteristic shape and predict that 



in a proton, the it-quark with helicity along the proton 
helicity carries about 2/3 of its longitudinal momentum 
(see Fig. || ). 

For the computation of the hard scattering amplitude 
Th (black circle in Fig. ||), the leading order PQCD contri- 
bution corresponds to the exchange of the minimum num- 
ber of gluons (two in the present case) between the three 
quarks. The number of diagrams grows rapidly with the 
number of elementary particles involved in the reaction 
(42 diagrams for the nucleon form factor, 336 diagrams in 
the case of real or virtual Compton scattering). Despite 
the large number of diagrams, the calculation of Th is 



a parameter free calculation once the scale Aqcd ~ 200 
McV in the strong coupling a s (Q 2 ) is given. Note that 
configurations with more than three valence quarks are a 
priori allowed but since this implies the exchange of more 
hard gluons, the corresponding contribution is suppressed 
by powers of 1/t. 




Fig. 6. Model distribution amplitudes for the nucleon : KS 
(upper figure) and COZ (lower figure), as function of the va- 
lence quark momentum fractions xi and X2 (xi + X2 + xz = 
1). 



There are two characteristic features of the Brodsky- 
Lepage model which are almost dir ect con sequences of 



QCD: the dimensional counting rules Bro73 and the con- 



servation of hadronic helicities | Bro81| . The latter feature 
implies that any helicity flip amplitude is zero and, hence, 
any single spin asymmetry too. The helicity sum rule is a 
consequence of utilizing the collinear approximation and 
of dealing with (almost) massless quarks which conserve 
their helicities when interacting with gluons. Whereas the 
dimensional counting rules are in reasonable agreement 
with experiment, the helicity sum rule seems to be vio- 
lated even at moderately large momentum transfers. The 
prevailing opinion is that these phenomena cannot be ex- 
plain ed in terms of perturbative QCD (see, for example, 
Ref. flSiv89|l ), but rather are generated by an interplay of 
perturbative and non-perturbative physics. 

An interesting aspect of real and virtual Compton scat- 
tering is that these are the simplest processes in which the 
integrals over the longitudinal momentum fractions yield 
imaginary parts. The reason is that, as in any scatter- 
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ing process, there are kinematical regions where internal 
quarks and gluons can go on their mass shell. The ap- 
pearance of imaginary parts to leading order in a s is a 
non-trivial prediction of PQCD, which should be tested 
experimentally. As discussed in [Gui98|, the (e, e'j) reac- 



tion with polarized incoming electrons seems to be a good 
candidate for this investigation. 

In contrast to the PQCD (or hard scattering ) approach 
to RCS, it was argued in Refs. Rad98a , Die99 that wide 
angle Compton scattering at accessible energies is described 
by a competing mechanism, in which the large momen- 
tum transfer is absorbed on a single quark and shared 
by the overlap of high-momentum components in the soft 
wave function. This so-called soft-overlap mechanism gives 
a purely real amplitude, therefore displaying a different 
signature than the PQCD amplitude. The transition from 
such a soft-overlap mechanism to the perturbative, hard 
scattering approach when increasing the momentum trans- 
fer is an open question for a reaction such as wide angle 
Compton scattering. It is hoped that future experiments 
can shed light on this transition. 



4.3 Results for RCS in PQCD 

The leading order PQCD prediction for RCS at large mo- 
mentum transfer has been calculated several times in the 



velop singularities : in the present case this number is 0, 1, 
2 or 3. Besides the trivial case of zero singularities which 
can be integrated immediately, the diagrams with one or 
two propagator singularities can be integrated by perform- 
ing a contour integration in the complex plane for one of 
the four integrations. For the most difficult case of three 
propagator singularities, it was found to be possible to 
evaluate the diagram by performing two contour integra- 
tions in the complex plane. In doing so, one achieves quite 
a fast convergence because the integrations along the real 
axis are replaced by integrations along semi-circles in the 
complex plane which are far from the propagator poles. 
This method was compared with the principal value in- 
tegration method, and the same r esults we re found up to 
0.1% for each type of singularity [ Vdh97b |. The principal 
value method was however found to converge much slower 
and is more complicated to implement, especially for the 
case with three singularities due to the coupled nature of 
the three principal value integrals. 



literature [|Far9q , |Far9li|Kro9l| , |Vdh97b| , |Broo00 |. 

The first step in such a calculation consists of evaluat- 
ing the 336 diagrams entering the hard scattering ampli- 
tude Tji for RCS. Next, the four-fold convolution integral 
of Eq. ( |37| ) has to be performed to obtain the Compton he- 
licity amplitudes. The numerical integration requires some 
care because the quark and/or gluon propagators can go 
on-shell which leads to (integrable) singularities. The dif- 
ferent numerical implementations of these singularities are 
probably the reason of the different results obtained in the 
literature. 



In Refs. [ Far90 , Far91 1 , the propagator singularities were 
integrated by taking a finite value for the imaginary part 
+ie of the propagator. The behavior of the result was then 
studied by decreasing the value of e. To obtain convergence 
with a practical numbe r of samples i n the Monte Carlo in- 
tegration performed in [ Far9C , Far9l|, the smallest feasible 
value for e was e w 0.005. In Ref. pCro91jl , the propagator 
singularities were integrated by decomposing the propa- 
gators into a principal value (off-shell) part and an on- 
shell pa rt. Both m ethods were implemented and compared 
in Ref. Vdh97b , and it was found that the +ie method 



yields differences of the order of 10% for every diagram as 
compared with the result of the principal value method. 
It is not surprising that, when summing hundreds of di- 
agrams, an error of 10% on every diagram can easily be 
amplified due to the interference between the diagrams. 

To have confidence in the evaluation of the convolu- 
tion of Eq. (|37j), the principal value integration method 
was compared in Ref. [ Vdh97b] with a third independent 
method. This third method starts from the observation 
that the diagrams can be classified into four categories 
depending upon the number of propagators which can de- 



Comparing the results of Ref. [ Vdh97b| with those 
of Ref. [Kro91|, a rather good agreement was found for 
all helicity amplitudes, except for the hclicity amplitude 
where both photon and proton helicities are positive, in 
which case both calculations differ strongly. Very recently, 
the PQCD calculation for RCS at large momentu m trans- 
fer has been recalculated again in Ref. |Broo00|, by also 
performing convolution integrals through contour integra- 



tions in the complex plane. The authors of Ref. |Broo00 



also find a strong difference with the results of [ Kro91| for 
the same helicity amplitude, where both photon and pro- 
ton helicities are positive. Furthermore, in the angular re- 
gion around 90°, where the PQCD fo rmalism is supposed 
to be applicable, the authors of Ref. [ BrooOOH fin d a good 
agreement with the calculations of Ref. |Vdh97b |, keeping 
in mind that there is an overall normalization uncertainty 
in these PQCD calculations for RCS, associated with a s 
and the valence quark wave function normaliza tion. The 
remaining differen ce betw een the results of Refs. [ Vdh97t ] 
and those of Ref. [ BrooOC ] seems to be isolated to a single 
hclicity amplitudes and appears for backward scattering 
angles. We therefore limit ourselves in the following dis- 
cussion to the results in the a ngular r egion around 90° 
where the calculations of Refs. |Vdh97t] and [BrooOC] are 
in good agreement, and which is the most relevant region 
for the PQCD calculation as it corresponds to the largest 
momentum transfer for a given photon energy. 

In Figs. and |, the PQCD calculation s for RC S are 
shown f or sev eral mo del DA's denoted as CZ [Che84], COZ 
|Che89| , KS jKin87| , and the asymptotic DA. 

The highest energy data which exist for real Comp- 
ton scattering were taken around 5 GeV and are shown in 
Fig. [?]. Although the energy at which these experiments 
were performed is probably too low to justify a PQCD 
calculation, the comparison with these data is neverthe- 
less shown in Fig. |7J for illustrative purposes. The nor- 
malization of the calculations shown at these very low 
scales corresponds to using a frozen coupling constant, 
with a s ss 0.5. 
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Fig. 7. Unpolarized Compton cross section on the proton for 
different nucleon DA's : KS (full line), COZ (dashed-dotted 
line), CZ (dashed line) and asym ptotic DA (dotted line). Cal- 
culations are from Ref. [Vdh97b|, where also the references of 
the data can be found. 



One first notices that the hard scattering amplitude 
for RCS has the s-dependence (T ~ s 



2 ) which leads 
da -- 6 for 



that is % r, s 



to the QCD scaling laws [Bro73| 

Compton scattering or VCS. The unpolarized real Comp- 
ton differential cross section (multiplied by the scaling fac- 
tor s 6 ) is shown in Fig. as function of the photon cm. 
angle. It is observed that the result with the asymptotic 
DA (~ 120x1X2^3) is more than one decade below the 
results obtained with the amplitudes KS, COZ, and CZ, 
motivated by QCD sum rules. The results with KS, COZ 
and CZ show a similar characteristic angular dependence 
which is asymmetric around 90°. Note that in the forward 
and backward directions, which are dominated by diffrac- 
tive mechanisms, a PQCD calculation is not reliable. Com- 
paring the results obtained with KS, COZ and CZ, one no- 
tices that although these DA's have nearly the same lowest 
moments, they lead to differences of a factor of two in the 
Compton scattering cross section. Consequently, this ob- 
servable is sensitive enough to distinguish between various 
distribution amplitudes, provided, of course, one is in the 
regime where the hard scattering mechanism dominates. 

In Fig. ||, the polarized Compton cross sections are 
shown for the two helicity states of the photon and for a 
target proton with positive helicity. One remarks that for 
all DA's there is a marked difference both in magnitude 
and angular dependence between the cross sections for the 
two photon helicities. Consequently, the resulting photon 
asymmetry S, defined as 



&(T,A = 1)- 



#(T,A = -1) 



where A is the helicity of the incoming photon and f de- 
notes a positive hadron helicity, changes sign for the DA's 
KS, COZ, and CZ for different values of <9 c . m . as shown 
in Fig. ||. It is seen that the asymptotic DA on the other 
hand yields a very large, negative a symmetr y around 90°. 
Therefore, it was suggested in Ref. [ Vdh97b] that the pho- 
ton asymmetry might be a particularly useful observable 
to distinguish between nucleon distribution amplitudes. 
The predicted sensitivity of the asymmetry to the nucleon 
DA can be used in the extraction of a DA fro m Compt on 
scattering data in the scaling region. In Ref. [Vdh97b], a 
procedure was outlined to extract a DA from Compton 
data in a model independent way by first expanding the 
DA in a set of basis functions and then using the angular 
information of the cross sections to fit the expansion coef- 
ficients. It was seen that the precision for these coefficients 
is greatly improved when one measures both unpolarized 
cross sections and photon asymmetries. 

7U + Pt ~> 7 + P (PQCD) 
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Fig. 8. PQCD calculations of the photon asymmetry on a 
polar ized prot on target for Compton scattering. Results (from 
Ref. ]ydh97b[ ) are shown for different DA's as indicated on 
the curves. 



A first dedicated experiment to measure the RCS dif- 
ferential cross section and the asymmetry of Eq. ( filj) for 
®c.m. around 90°, an d for a real photon energy of 6 GeV, 
is planned at JLab Woj99| . In particular, it will be in- 
teresting to see if one approaches the PQCD result at 
these "lower" energies, and to study the interplay with 
soft-o verlap type con tributions for RCS as proposed in 
Refs. ]Rad98a| , pic99fl . 

RCS" experiments using a real photon energy in the 15 
GeV range, might be f easible e.g. at t he HERA ring in 
the foreseeable future | d'Ho96 , D'An97 | and might open 
up prospects to study the nucleon valence wave function 
in a direct way. 
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5 Deeply virtual Compton scattering and 
skewed parton distributions 

5.1 Introduction 

Much of the internal structure of the nucleon has been 
revealed during the last two decades through the inclu- 
sive scattering of high energy leptons on the nucleon in 
the Bjorken -or "Deep Inelastic Scattering" (DIS)- regime 
(where the photon virtuality Q 2 is very large, and xb — 
Q 2 /2p.q finite) . Unpolarized DIS experiments have mapped 
out the quark and gluon distributions in the nucleon, while 
polarized DIS experiments have shown that only a small 
fraction of the nucleon spin is carried by the quarks. This 
has stimulated new investigations to understand the nu- 
cleon spin. 

With the advent of the new generation of high-energy, 
high-luminosity lepton accelerators combined with large 
acceptance spectrometers, a wide variety of exclusive pro- 
cesses in the Bjorken regime are considered as experi- 
mentally accessible. In recent years, a unified theoreti- 
cal description of such processes has emerged through 
a formalism introducing a new type of parton distribu- 
tions, co mmonly deno ted as skewed parton distributions 
(SPD's) |lJi97t|Rad96a|| . These SPD's are generalizations 
of the parton distributions measured in DIS. It has been 
shown that these SPD's, which parametrize the struc- 
ture of the nucleon, allow one to describe, in leading or- 
der perturbative QCD (PQCD), various exclusive pro- 
cesses in the near forward direction, where the momentum 
transfer to the nucleon is small. Such non-forward pro- 
cesses were alrea dy cons i dered i n the literature a l onger 
time ago, see e.g. JWat82| , pa7^ , |Dit88| , [lai93||Mul94| . The 
most promising of these non-forward hard exclusive pro- 
cesses are deeply virtual Compton scattering (DVCS) and 
longitudinal clcctroproduction of vector or pseudoscalar 
mesons at large Q 2 . 



5.2 Definitions and modelizations of skewed parton 
distributions 

The leading order PQCD diagrams for DVCS and hard 
meson electroproduction are of the type as shown in Fig. [| 
The hard scale in Fig. [9] is the photon virtuality Q 2 , which 
should be large (of the order of several GeV 2 ), so as to be 



in the Bjorken regime. It has been proven | Ji97, Rad96a| 
that the leading order DVCS amplitude in the forward di- 
rection can be factorized in a hard scattering part (which 
is exactly calculable in PQCD) and a soft, nonpertur- 
bative nucleon structure part as illustrated on the left 
panel of Fig. ^. The nucleon structure information can 
be parametrized, at leading order, in terms of four (quark 
helicity conserving) generalized structure functions. These 
functions are the SPD's denoted by H, H, E, E which de- 
pend upon three variables : x, £ and f. The light-cone 
momentum ^| fraction x is defined by k + = xP + , where k 

3 using the definition a = 1/a/2(i° ± a 3 ) for the light-cone 
components 





A A 



Fig. 9. Leading order diagrams for DVCS (left) and for longi- 
tudinal electroproduction of mesons (right). 



is the quark loop momentum and P is the average nucleon 
momentum (P = (p+p')/2, where p{p') are the initial (fi- 
nal) nucleon four- momenta respectively) . The skewedness 
variable £ is defined by A + — — 2£ P + , where A = p' — p is 
the overall momentum transfer in the process, and where 
2£ — > Xb/(1 ~ xb/2) in the Bjorken limit. Furthermore, 
the third variable entering the SPD's is given by the Man- 
delstam invariant t = A 2 , being the total squared momen- 
tum transfer to the nucleon. In a frame where the virtual 
photon momentum q^ and the average nucleon momen- 
tum P M are collinear along the z-axis and in opposite di- 
rection, one can parametrize the non-perturbative object 
in the lower blobs of Fig. || as : 



P+ 
27 



dy e 



ixP + y 



(p \M-y/2)i> a (y/2)\p) 



\{(Tf-)c0 [H q (x,£,t) N(p') 7 + N(p) 



y + =y±=o 



+ E q (x,^t) N(p )ia 



+K 



A„ 



N(p) 



-(757" 



)uf3 



2tojv 

H q {x,i,t) N{p)^ l5 N{p) 

+ EO(x,Z,t) N(p') l5 ^—N(p) 
2m,N 



(42) 



where V is the quark field, TV the nucleon spinor and 
totv the nucleon mass. The Ihs of Eq. ( j42] ) can be in- 
terpreted as a Fourier integral along the light-cone dis- 
tance y~ of a quark-quark correlation function, represent- 
ing the process where a quark is taken out of the ini- 
tial nucleon (with momentum p) at the space-time point 
y/2, and is put back in the final nucleon (with momen- 
tum p') at the space-time point —y/2. This process takes 
place at equal light-cone time (y + — 0) and at zero trans- 
verse separation (yj_ = 0) between the quarks. The re- 
sulting one-dimensional Fourier integral along the light- 
cone distance y~ is with respect to the quark light-cone 
momentum xP + . The rhs of Eq. (|42|) parametrizes this 
non-perturbative object in terms of four SPD's, according 
to whether they correspond to a vector operator (7" ) aj g 
or an axial-vector operator (757~) a ,3 at the quark level. 
The vector operator corresponds at the nucleon side to a 
vector transition (parametrized by the function H q , for a 
quark of flavor q) and a tensor transition (parametrized by 



16 



Marc Vanderhaeghen: Real and Virtual Compton Scattering off the Nucleon 



the function E q ). The axial- vector operator corresponds 
at the nucleon side to an axial-vector transition (function 
H q ) and a pseudoscalar transition (function E q ). 

In Fig. O, the variable x runs from -1 to 1. Therefore, 
the momentum fractions [x + £ or x — £ ) of the active 
quarks can either be positive or negative. Since positive 
(negative) momentum fractions correspo nd to quarks (an- 
tiquarks), it has been noted in [Rad96a| that in this way, 

when x > £ 
-£ both 



one can identify two regions for the SPD's 
both partons represent quarks, whereas for x < 
partons represent antiquarks. In these regions, the SPD's 
are the generalizations of the usual parton distributions 
from DIS. Actually, in the forward direction, the SPD's H 
and H reduce to the quark density distribution q(x) and 
quark helicity distribution Aq{x) respectively, obtained 
from DIS : 



H q {x,0,0) = q(x) , H q (x,0,0) = Aq(x) 



(43) 



The functions E and E are not measurable through DIS 
because the associated tensors in Eq. ( p2[) vanish in the 
forward limit (A — ► 0). Therefore, E and E are new lead- 
ing twist functions, which are accessible through the hard 
exclusive electroproduction reactions, discussed in the fol- 
lowing. 

In the region — £ < x < £, one parton connected to 
the lower blob in Fig. ^ represents a quark and the other 
one an antiquark. In this region, the SPD's behave like a 
meson distribution amplitude and contain completely new 
information about nucleon structure, because the region 
— £ < x < £ is absent in DIS, which corresponds to the 
limit £ -> 0. 

Besides coinciding with the quark distributions at van- 
ishing momentum transfer, the skewed parton distribu- 
tions have interesting links with other nucleon structure 
quantities. The first moments of the SPD's are related to 
the elastic form factors (FF) of the nucleon through model 
independent sum rules. By integrating Eq. ( [l2| ) over x, one 
obtains the following relations for one quark flavor : 



+i 



dxH q (x,^t) 



dxE q (x,£,t) 



l 

+ i 



dxH q (x,^t) 



dxE q (x,£,t) 



9 Q A (t), 
h\{t). 



(44) 



The elastic FF for one quark flavor on the rhs of Eqs. ( |44|) 
are related to the physical ones (restricting oneself to u, d 
and s quark flavors) as : 



= 2Ff 



F™ + Fl , Ff = 2F{ 1 



Ff 



F?, (45) 



where Ff and F" are the usual proton and neutron Dirac 
FF respectively, and where Ff is the strangeness form fac- 
tor. Relations similar to Eq. ( fr5|) hold for the Pauli FF F| . 



For the axial vector FF one uses the isospin decomposi- 
tion : 



1 



1 



9a = 2 9A + 2 9a 



1 



1 



--j9A + ^ 9 A 



(46) 



where 3^(9^4) ar e the isovector (isoscalar) axial FF respec- 
tively. Similar relations exist for Ha- The isovector axial 
form factor g& is known from experiment, with (?a(0) ~ 
1.267. The induced pseudoscalar form factor Ha contains 
an important pion pole contribution, through the partial 
conservation of the axial current (PCAC). 

A lot of the recent interest and activi ty in this field 
has been triggered by the observation of | Ji97 | that the 
SPD's may shed a new light on the "spin-puzzle". Start- 
ing from a (color) gauge-invariant decomposition of the 



nucleon spin : 1/2 = J q + J g 



where J q and J g are the 



total quark an d glu on angular momentum respectively, it 
was shown in | Ji97 | that the second moment of the unpo- 
larized SPD's at t = gives 



dxx [H q {x^,t = 0) + E q {x,£,t = 0)}, (47) 



and this relation is independent of £. The quark angular 
momentum J q decomposes as : J q — AS/ 2 + L q , where 
AU/2 and L q are the quark spin and orbital angular mo- 
mentum respectively. As AS is measured through polar- 
ized DIS experiments, a measurement of the sum rule of 
Eq. ( |47| ) in terms of the SPD's, provides a model indepen- 
dent way to determine the quark orbital contribution L q 
to the nucleon spin. 

Ultimately, one wants to extract the SPD's from data, 
but in order to evaluate the electroproduction observables, 
and to study their sensitivity to the new ph ysics, o ne 
needs an educated guess for the SPD's. In Ref. | Vdh99 |, a 
model for the SPD's was constructed using a ^-dependent 
prod uct ansat z (for the double distributions introduced in 
Ref. |Rad98b|) of a quark distribution and an asympt otic 
"meson-like" distribution amplitude (see Ref. [ Vdh9£ ] for 
more details). F or the q uark distributions, the MRST98 
parametrization | Mar98| is used as input. The ^-dependence 
of the model for the SPD's is given by the corresponding 
FF (Dirac form factor for H, axial form factor for H), 
so that the first moments of the SPD's are satisfied by 
construction. As an example, the (i-quark SPD (formerly 
also denoted as off- forward parton distribution (OFPD)), 
using the above described ansatz, is shown in Fig. [l^. 

One observes from Fig. [l^ the transition from a quark 
distribution (£ = 0) to a meson distribution amplitude 
(£ = 1). Model calculation of the SPD's are currently pos- 
sible within the QCD chir al mode ls for intermediate xb- 
In particular, a calcul ation [ | Pct98 in the chiral quark soli- 
ton model (see Ref. | Chr96| for a review) found a strong 
dependence of the SPD's on £ and fast "crossovers" at 

= £. Such behavior is related to the fact that the SPD's 
in the region — £ < x < £ have properties of meson dis- 
tribution amplitudes. In particular for the SPD H, this 
can be seen as being due to a scalar-isoscalar two-pion ex- 
change contribution [Pol9£], indicating that the SPD's are 
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Fig. 10. £ dependence of the SP D H d at t = with the ansatz 
(bas ed on th e MRST98 [Mar9£] quark distributions) used in 
Ref. | Vdh99 . Upper panel : valence d-quark SPD, lower panel : 
total d-quark SPD. The thin lines (£ = 0) correspond with the 
ordinary d-quark distributions. 



qualitatively a richer source of nucleon structure informa- 
tion than ordinary parton distributions. One may expect 
that eventually it will be possible to calculate SPD's for 
intermediate xb using lattice QCD. 



\ / dx 



x — £ + ie x + £, — it 



H p DVCS (x,U) N(p'h.nN(p) 
+ E p DVCS {x, £, t) N(p) t a KX ^Nun 



+i 



dx 



1 



x — I; + ie x + I; — ie 



_ , A ■ n 
+ E p DVCS (x,^t)N(p ) 75 - N(p) 



2m 



N 



(49) 



On the rhs of the DVCS tensor of Eq. (|49|), the SPD's 
H, H, E, E enter in a convolution integral over the quark 
momentum fraction x. This is a qualitative difference com- 
pared with the DIS amplitude, where one is only sensitive 
(through the optical theorem) to the imaginary part of 
the for ward do uble virtual Compton amplitude. We refer 
to Ref. Gui98 for details and for the formalism to calcu- 
late DVCS observables starting from the DVCS tensor of 
Eq. ©. 

The leading order DVCS amplitude corresponding to 
Eq. (|49|), is exactly gauge invariant with respect to the 
virtual photon, i.e. q v H£ V q dvcs = 0. However, electro- 
magnetic gauge invariance is violated by the real photon 
except in the forward direction. This violation of gauge in- 
variance is a higher twist (twist-3) effect compared to the 



Since q H^ u 



leading order term H^ Q DVCS . h^^l.o.dvcs 
A±, an improved DVCS amplitude linear in A± has been 
proposed in Ref. [ Gui98| ] to restore gauge invariance (in 
the nonforward direction) in a heuristic way : 



5.3 Leading order amplitudes and observables for 
DVCS and hard meson electroproduction 

The leading order (DO.) DVCS amplitude in the forward 
direction is given ]Ji9 j[ by the handbag diagram shown on 
the left panel of Fig. |9] (the crossed diagram which is not 
shown is also understood). A for mal factoriza tion proof 
for DVCS has been given in Rcfs. [jl98a| , |Col99C . 

To calculate the DVCS amplitude in the Bjorken regime, 
it is natural to express the momenta in the process (q^ of 
the virtual photon, q ,fl of the real photon, and P M denoting 
the average nucleon momentum) in terms of the lightlike 
vectors 



2^ = ^(1,0,0,1) 



1 



P+V2 



(1,0,0,-1). (48) 



Using the parametrization of Eq. ( ^2[) for the bilocal quark 
op erator the L.O. DVCS tensor H^ DVC s (defined e.g. 
in |Gui98 |) follows from the handbag diagrams as : 



n L.O. DVCS 



^DVCS ~ n L.O.DVCS 



(p ■ q 



L.O. DVCS 



(50) 

DVCS amplitude 



leading to a correction term to the L.O. 
of order 0(A ± /Q). 

Very recently, the gauge invariance of the DVCS ampli- 
tude was addressed in much more detail in several works 



[|Ani0C| , |Pen00| , |Bel00b| . It was found that the twist-3 terms 
which restore gauge invariance (to twist-4 accuracy) in- 
volve two contributions. First there are terms proportional 
to the twist-2 SPD's of Eq. (|9|), which were found to 
completely coincide with the improved DVCS amplitude 
of Eq. (J50j). In addition, there are terms whi ch are charac- 
terized by new 'transverse' SPD's (see Refs. AniOO] , PenOO 



p3cl00b| for details). These latter functions are suppressed 
by one power 1/Q compared with the contribution of the 
twist-2 SPD's in DVCS cross sections, and could in prin- 
ciple be separated by measuring DVC S obse rvables over 
a sufficiently large Q 2 range (see e.g. [Die97| for tests of 
the handbag approximation to DVCS). In view of current 
DVCS experiments which are performed or planned at Q 2 
in the few GeV 2 range only, the numerical importance of 
those additional contributions remains to be investigated. 
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Besides the DVCS process, a factorization proof was 
also given for the L.O. mes on electroprodu ction ampli- 
tudes in the Bjorken regime [Col97,|Rad96b|, which is il- 



lustrated on the right panel of Fig. |9j. This factorization 
theorem only applies when the virtual photon is longitudi- 
nally polarized. In the valence region, the L.O. amplitude 
A4 L for meson production by a longitudinal photon con- 
sists of evaluating Fig. || (right panel, where only one of 
the four L.O. diagrams is shown) with the one-gluon ex- 
change diagrams as hard scattering kernel. In this way, 
the L.O. expressions for p Q L (longitudinally polarized vec- 
tor mes on) and ir° ele ctropro duction were calculated in 
[ Vdh9SH (see also Ref. flMan98| ) as : 



factors. For DVCS on the proton, the combination is 



= -H 1 



-H a 



-H s 



(53) 



and similarly for H, E and E. For electroproduction of 
p° and 7T° on the proton, the isospin structure yields the 
combination 



o 

-H d 

3 



(54) 



(55) 



M 



PL 



Mi 



4 1 
" 6 9Q 



dz 



1 
2 

(Ana, 



dx 



- £ + it x + £ — it 



+E p (x,Z,t)N(p)ia K 



2lTlN 



N(p) ,(51) 



4 1 



dx 



dz 



i 



i 



x — £ + it x + £ — ? 
(A<ira s ) {i% (x,£,t)N(p' )l-n l5 N{p) 

+E%(x,tt)N(p)j 5 ^N(p) 



(52) 



and similar for E and E. Corresponding relations for the 



7r and 7 7 chann els can be found in Refs. [Fra99 



Man99a , Man99t , Vdh99 1 . Therefore, the measurements of 
the different meson production channels are sensitive to 
different combinations of the same universal SPD's, and 
allow us to perform a flavor separation of the SPD's, pro- 
vided one is able to deconvolute the SPD's from the lead- 
ing order amplitudes. 

Some representative results for DVCS and meson elec- 
troproduction observables using the ^-dependent ansatz 
for the SPD's, are shown in t he following. More de tailed 
results can be found in Refs. ]Vdh98| , |Gm98| , |Vdh99| | . 

Before considering the extraction of the SPD's from 
electroproduction data, it is compulsory to demonstrate 
that the scaling regime has been reached. In Fig. |ll|, the 
forward longitudinal electroproduction cross sections are 
shown as a function of Q 2 and the L.O. predictions are 
compared for different mesons. The L.O. amplitude for 
longitudinal electroproduction of mesons was seen to be- 
have as 1 /Q, leading to a 1 /Q 6 scaling behavior for da^ / dt. 



where a s is the QCD coupling constant. Because the quark 
helicity is conserved in the hard scattering process, one 
finds the interesting result that the vector meson elec- 
troproduction amplitude depends only on the unpolar- 
ized SPD's H and E, whereas the pseudoscalar meson 
electroproduction amplitudes depend only on the polar- 
ized SPD's H and E. In contrast, the DVCS amplitude 
of Eq. ( f49|) depends on both the unpolarized and polar- 
ized SPD's. Another difference from DVCS, is the fact 
that the meson electroproduction amplitudes require ad- 
ditional non-perturbative input from the meson distribu- 
tion amplitudes ^p(z) and ^(z) respectively, for which 
the as ympto tic forms are taken in the calculations. From 
Eqs. (51,52), one furthermore sees that the L.O. 



longi- 
tudinal amplitudes for meson electroproduction behave as 
l/Q- At large Q 2 , fixed xb and fixed t, this leads to a 1/Q 6 
behavior for the longitudinal cross section daL/dt, which 
provides an experimental signature (scaling) of the lea ding 
order mechanism. Expressions analogous to Eqs. (plj ^2[) 
have also been worked out for the charged meson chan- 
nels p ^- , 7r^ as well as for the uj, <j> a nd r\ channels (see 



Refs. [ Fra9£ , Man99a , Man99b , Vdh99 1 for details) 



> 
o 



II 10 2 



— 




According to the considered reaction, the SPD's enter 
in different combinations due to the charges and isospin 



10 

Q 2 (GeV 2 ) 

Fig. 11. Scaling behavior of the L.O. predictions for the for- 
ward differential electroproduction cross section on the pro- 
ton, for vector mesons (left panel) and ps eudoscalar mesons 
(right panel), as calculated in Ref. [Vdh99]. For the tt + chan- 
nel, the pion pole contribution (full line, (7r + ) po ; e ) is shown 
separately from the H contribution (dashed line, (7r + )„. p .). 
Also shown is the scaling behavior of the forward transverse 
cross section dar/dt for the leading order DVCS cross section 
(dashed-dotted line in left panel). 
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By comparing the different vector meson channels in 
Fig. O, one sees that the p\ channel yields the largest 
cross section. The col channel in the valence region (xb ~ 
0.3) is about a factor of 5 smaller than the p\ channel, 
which is to be compared with the ratio at small xb (in 
the diffractive regime) where p° : ui = 9 : 1. The p\ chan- 
nel, which is sensitive to the isovector combination of the 
unpolarized SPD's, yields a cross section comparable to 
the ojl channel. The p\ channel is interesting as there 
is no competing diffractive contribution, and therefore al- 
lows to test directly the quark SPD's. The three vector 
meson channels pj, uol) are highly complementary 
in order to perform a flavor separation of the unpolarized 
SPD's H u and H d . 

A dedicated ex perimen t is planned at JLab at 6 GeV 
in the near future jGuid9£ ] to investigate the onset of the 



have been proposed for DVCS ]d'Ho99a| and meson elec- 
troproduction [Poc9£|. 



scaling behavior for p v L electroproduction in the valence 
region (Q 2 » 3.5 GeV 2 , x B ~ 0.3). 

For the pseudoscalar mesons which involve the polar- 
ized SPD's, one remarks in Fig. [ll] the prominent contri- 
bution of the charged pion pole to the ir + cross section. 
For the contribution proportional to the SPD H, it is also 
seen that the ir° channel is about a factor of 5 below the 
7T + channel due to isospin factors. In the 7r° channel, the 
u- and d-quark polarized SPD's enter with the same sign, 
whereas in the 7r + channel, they enter with opposite signs. 
As the polarized SPD's are constructed here from the cor- 
responding polarized parton distributions, the difference 
between the predictions for the tt° and ir + channels re- 
sults from the fact that the polarized d-quark distribution 
is opposite in sign to the polarized it-quark distribution. 
For the ij channel, the ansatz for the SPD H based on the 
polarized quark distributions yields a prediction compa- 
rable to the 7T° cross section. 

For the 7 leptoproduction in the few GeV beam en- 
ergy range, the cross section is domin ated by the Bethe- 
Hcitler (BH) proce ss (see Ref. [ Gui98| ). However, it was 
suggested in Ref. | Gui98| that an exploration of DVCS 
might be possible if the beam is polarized. The electron 
single spin asymmetry (SSA) does not vanish out of plane 
due to the interference between the purely real BH process 
and the imaginary part of the DVCS amplitude. Because 
the SSA measures the imaginary part of the DVCS ampli- 
tude, it is directly proportional to a linear combination of 
the SPD's along the line x = £. In fact, the SSA maps out 
an 'envelope' function, e.g. H(x — £,£,i), as shown e.g. 
in Fig. [12] for the valence down quark SPD in the ansatz 
corresponding with Fig. [l^. 

In Fig. O, it is shown that the SSA yields a sizeable 
asymmetry for JLab kinematics, and displays a sensitivity 
to the ^-dependent shape of the SPD's. An experiment to 
measure the SSA for DV CS has very recently been pro- 
posed at JLab at 6 GeV (Bert9|. The SSA fo r DVCS is 
at present also measured at HERMES [ Arm99 ]. 

Going up in energy, the increasing virtual photon flux 
factor boosts the DVCS part of the 7 leptoproduction 
cross section, making it more important compared to the 
BH contribution. This provides a nice opportunity for 
COMPASS at 200 GeV beam energy, where experiments 
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Fig. 12. Valence d-quark contribution to the SPD H d at t = 
for different values of £ as indicated on the curves, calculated 
with the ansatz as in Fig. [To| The thin (lower) solid curve (£ = 
0) corresponds to the ordinary d-quark distributions, whereas 
the thick (upper) solid curve to the envelope function Hy(x = 
£, £, t = 0) as measured through the SSA. 
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Fig. 13. Single spin asymmetry for DVCS in JLab kinematics. 
A comparison is shown of the a symme tries for a ^-independent 
ansatz for the SPD's as in Ref. |]Gui98 1 (dashed-dotted curves), 
and a ^-dependent ansatz for the SPD's as in Fig. (full 
curve) . 
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5.4 Extension to hard exclusive electroproduction of 
decuplet baryons 

In the previous section, the main focus were the reactions 
7* + N -> 7 + N' and jl + N -> M + N' with M a meson, 
and where AT' is an octet baryon. One of the intriguing 
questions of medium-energy QCD dynamics is the differ- 
ences and similarities in the structure of baryons belonging 
to the different SU(3)f multiplets. In particular, a naive 
constituent quark model predicts that they are similar, 
while there are suggestions that due to a strong attrac- 
tion between the quarks in the spin-isospin zero channel, 
diquark correlations should be impor tant in the baryon 
octet but not in the decuplet [ 3ch98| . At the same time 
the chiral models suggest that in the limit of a large num- 
ber of colors (large N c ) of QCD, which is known to be a 
useful guideline, nucleons and A isobar s are different ro- 
tational excitations of the same soliton [ Adk83, Dia88|. 

For these studies, the potential of the process 72+A" — » 
7r + A as well as th e DVC S process 7* + N — > 7 + A, 
was explored in Ref. |FraOC ] . In addition, the study of the 
processes with production of decuplet baryons has also a 
practical usefulness, because in the experiments with low 
resolution in the mass of the recoiling system (AM « 300 
MeV for HERMES in the current set-up), the estimates of 
A production are necessary to extract the N —> N SPD's 
from such data. 

In Ref. [ FraOOf , a new set of SPD's were introduced for 

/o\ 

the axial N —> A (isovector) transition, denoted as C\ , 
which enter into it A electroproduction : 



.7 p- 



oi ■ - : al ,W7t " « 1 : 0.5 : 1.25 : 0.8. 

Besides the cross section <7l, the second observable in- 
volving only longitudinal amplitudes and being a leading 
order observable for hard exclusive meson electroproduc- 
tion, is the single spin asymmetry, for a proton target po- 
larized perpendicular to the reaction plane (or the equiv- 
alent recoil polarization observable) Fra99| . These trans- 



>7T p 



P+_ 

2tt 



dy-e ixP+ *- (A+ \$(-y/2)jry 5 \N) 



v+=y±=o 



Ap(n-A) 



N 



N(p), (56) 



where the same notations are used as before, and where 
ifiP(p') is the Rarita-Schwinger spinor for the A isobar. In 
Eq. ([56|) , the ellipses denote other contributions which are 
suppressed at large N c . For the A^ — * A DVCS process, 
besides the axial SPD's, also v ector SPD's enter which 
were also defined in Ref. [ Fra00| . 

The observation that in the large A" c limit, the nucleon 
and A are rotational excitations of the same classical soli- 
ton object, allows us to derive a number of relations be- 



tween N ^ N and AT -> A SPD's. For c} 3) and C. 



<3) 



these have the form [ Fra99| 



(57) 
(58) 



in terms of the (isovector) N ^ N SPD = H u - H d , 
and analogously for E^ 3 \ 

Using the large N c relations of Eq. (^), one can eas- 
ily derive the relations between the different cross sections 



for charged pion production as a L 



7*p — >7r^n _ 7*P — >tt^ A° 



vers e sp in asymmetries for n + n and tt + A 
Fig. |l 
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are shown in 
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Fig. 14. Transverse spin asymmetries for the longitudinal elec- 
troproduction of 7r + n and 7r + Z\°, a t different values of t (in 
(GeV/c) 2 ). Figure from Ref. §raO0|. 



: a 



It is obvious from Fig. [l4|, that large transverse spin 
asymmetries are predicted for these processes, related to 
the peculiar feature of chiral QCD. As a consequence, in- 
vestigations of these processes can provide unique tests of 
the soliton type approach to baryon structure. The spin 
asymmetries are likely to be less sensitive to higher twist 
effects and hence can be explored already using the HER- 
MES detector and JLab at higher energies. Furthermore, 
the study of these processes would allow one to make a 
more reliable separation of the n pole contribution in the 
electroproduction of pions, which is mandatory for the 
measurement of the pion elastic form factor at higher Q 2 . 



5.5 Power corrections to the leading order amplitudes 

When measuring hard electroproduction reactions in the 
region Q 2 w 1 — 20 GeV 2 , there arises the question of 
the importance of power corrections to the leading order 
amplitudes, i.e. how fast does one approach the scaling 
regime predicted by the L.O. amplitudes. One source of 
power corrections is evident from the structure of the ma- 
trix element of Eq. ( fl2|) which defines the SPD's, where 
the quarks are taken at zero transverse separation. This 
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amounts to neglect, at leading order, the quark's trans- 
verse momentum compared with its large longitudinal (+ 
component) momentum. A first estimate of these correc- 
tions due to the quark's in trinsic transverse momentum 
has been obtained in Ref. |Vdh99|, which is referred to 



for details. This correction is known to be important for a 
successful description at the lower Q 2 values of the 7r°7*7 
transition FF, for which data exist in the range Q 2 w 1—10 
GeV 2 . For the pion elastic FF in the transition region be- 
fore asymptotia is reached, the power corrections due to 
both the transverse momentum dependence and the soft 
overlap mechanism (i.e. the process which does not pro- 
ceed through one-gluon exchange) are quantitatively im- 
portant. The result for the pion elastic FF is shown in 
Fig. [if], where it is seen that the leading order PQCD re- 
sult is approached only at very large Q 2 . The correction 
including the transverse momentum dependence gives a 
substantial suppression at lower Q 2 (about a factor of two 
around Q 2 « 5 GeV 2 ). At these lower Q 2 values, the in- 
clusion of the transverse momentum dependence renders 
the PQCD calculation internally consistent. 

PION ELECTROMAGNETIC FORM FACTOR 
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O 
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Q ( GeV ) 



Fig. 15. Results for the ir electromagnetic form factor of the 
leading order PQCD prediction without (dotted line) and with 
(dashed-dotted line ) inclusion of the corrections due to intrin- 
sic transverse momentum dependence. The dashed curve shows 
the result for the soft overlap contribution and the total result 
(full line) is the sum of the dashed and dashed-dotted lines. 
Figure from Ref. [Vdh99], where the references of the data can 
also be found. 



These form factors were taken as a guide in Ref. [ Vdh99| 
to estimate the corrections due to the parton's intrinsic 
transverse momentum dependence in the DVCS and hard 
meson electroproduction amplitudes. 

Although experimental data for p® L electroproduction 
at larger Q 2 do not yet exist in the valence region (xb ~ 
0.3), the reaction j L p — > p° L p has been measured at 
smaller values of xb- Therefore, in Fig. [l(] the calcula- 
tions are compared to those data, in order to study how 



the valence region is approached, in which one is sensitive 
to the quark SPD's. For the purpose of this discussion, we 
call the mechanism proceeding through the quark SPD's 
(i.e. the right panel of Fig. [)]), the Quark Exchange Mech- 
anism (QEM). Besides the QEM, p° electroproduction at 
large Q 2 and small xb proceeds predominantly through a 
perturbative tw o-gluo n exchange mechanism (PTGEM) as 
studied in Ref. | Fra96 |. To compare to the data at interme- 
diate Q 2 , the power corrections due to the parton's intrin- 
sic transverse momentum dep endence were implemented 
in both mechanisms (see Ref. [ Vdh99| for details), which 
gives a significant reduction at the lower Q 2 . The results 
are compared with the data in Fig. |l6], showing that the 
PTGEM explains well the fast increase of the cross section 
at high cm. energy (W), but substantially underestimates 
the data at the lower energies. This is where the QEM is 
expected to contribute since xb is then in the valence re- 
gion. The results including the QEM describe the change 
of behavior of the data at lower W quite nicely. 

Y* + p^PL + p 




10 10 

W (GeV) 

Fig. 16. Longitudinal forward differential cross section for pi 
electroproduction. Calculations compare the quark exchange 
mechanism (dotted lines) with the two-gluon exchange mech- 
anism (dashed lines) and the sum of both (full lines). Both 
calculations include the corrections due to intrinsic transverse 
momentum dependence. The data are from NMC (triangles), 
E665 (solid circles), ZEUS 93 (open circles) and ZEUS 95 (open 
squares). Figure from Ref. [Vdh99 , where the references of the 
data can also be found. 

Recently, p\ data have been obtained by the HERMES 
Colla boratio n for Q 2 up to 5 GeV 2 and around W as 5 
GeV [ AirOO . These data show a clear dominance of the 
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QEM in the intermediate W range as predicted in | Vdh98 



Vdh99| . The model ansatz for the SPD's of Ref. [|Vdh99 | 



gives a fairly good agre ement with these longitudinal p u 
electroproduction data | AirO0|] . 



5.6 Perspectives and outlook 

In order to extract SPD's from forthcoming data, the Q 2 
evolution of the SPD's has alr eady b een worked out at 
the next-to-leading order level [Bel98|, which shows that 
the Q 2 evolution of the SPD's interpolates between the 
evolution of the parton distributions and the evolution of 
the distribution amplitudes. Also radiative corrections to 
the coefficient functio ns have been calculated recently in 
next-to-leading order [BelOOa]. 

A major open theoretical question in this field is how 
the SPD's can be deconvoluted from the leading order 
amplitudes. Suitable parametrizations of the SPD's, in- 
corporating all physical constraints, might be one avenue 
to tackle this problem. In absence of a solution to this 
problem, one has to resort to model calcuations or edu- 
cated guesses for the SPD's in order to compare with the 
data. 

On the experimental side, it is clear that new and accu- 
rate data are needed for various exclusive channels at large 
Q 2 in the valence region, where the quark exchange mecha- 
nism dominates. Several experimen ts are being performed 
or are planned or proposed at JLab [ Suid98 , Bert9£ , BertOC ] 
HERMES and COMPASS poc99Hd'Ho99af ■ Looking some- 
what further into the future, the measurement of hard ex- 
clusive reactions will be one of the cen tral the mes for the 
planned upgrade of JLab to 12 GeV |BurOC ] . A facility 
with high luminosity combined with an intermediate en- 
ergy of around 25 GeV, such as e.g. the ELFE project 



[Bur99 , will be a dedicated facility to measure these ex- 
clusive reactions at high momentum transfer and to map 
out these new SPD's in detail. Although such exclusive 
experiments at large Q 2 are quite demanding, the fun- 
damental interest of the SPD's justifies an effort towards 
their experimental determination. 



6 QED radiative corrections to virtual 
Compton scattering 

6.1 Introduction 

We have discussed in section || how VCS below pion pro- 
duction threshold, allows us to access generalized polariz- 
abilities of the proton. Furthermore, we have seen in sec- 
tion U that VCS in the Bjorken regime determines general- 
ized parton distributions of the nucleon. In both regimes, 
experiments are either being done or planned for the near 
future. In order to extract the nucleon structure informa- 
tion of interest from the ep — > ep7 reaction, especially in 
those kinematical situations where the Bethe-Heitler pro- 
cess is not negligible, it is indispensable to have a very 
good understanding of the QED radiative corrections to 



the ep — > ep7 reaction. 

The ep — > epr/ reaction is particular in comparison 
with other electron scattering reactions, because the pho- 
ton can be emitted from both the proton side (this is the 
VCS process which contains the nucleon structure infor- 
mation of interest) or from one of the electrons (which is 
the Bethe-Heitler process). The radiative corrections ob- 
tained from the Bethe-Heitler process differ formally from 
the case of elastic electron scattering. 



6.2 Results for the QED radiative corrections to VCS 

The first order QED radiative corre ctions to the ep — ► 
ep7 reaction were calculated in Ref. [VdhOOa]. The one- 



loop virtual radiative corrections have been evaluated by a 
combined analytical-numerical method. Several tests were 
performed to cross -check t he numerical method used. It 
was shown in Ref. [ VdhOOa how all IR divergences cancel 
when adding the soft-photon emission processes. Further- 
more, a fully numerical method was presented for the pho- 
ton emission processes where the photon energy is not very 
small compared with the electron energies, which makes 
up the radiative tail. 
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Fig. 17. Differential ep — > ep'y cross section for MAMI kine- 
matics at q = 111.5 MeV/c. Dashed-dotted curve : BH + 
Born contribution, dashed curve : BH + Born + virtual ra- 
diative correction, full curve : BH + Born + total radiative 



correction. Figure from Ref. [VdhOOa 



Fig. [I?] shows as representative result the effect of the 
radiative corrections on the VCS differential cross section 



for the MAMI VCS experiment [Roc00| at an outgoing 
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photon energy of q' = 111.5 MeV (we refer to |Vdh00a| 
for all details and more results). It is seen that the to- 
tal effect of the radiative corrections in these kinematics 
is a reduction of the BH+Born cross section of the order 
of 20%. The effect of the radiative corrections was also 
confirmed by the experimental results at very low energy 
of the outgoing photon (q' = 33 MeV), where the effect 
of the GP's is negligible. From the difference between the 
radiatively corrected data and the BH + Born result, the 
two values of Eq. (|2^) for the combinations of the pro- 
ton's GP's at Q 2 ~ 0.33 GeV 2 have been extracted in 

Ref. URocOOt 

In Rcf. Vdh00a|, calculations of the VCS radiative 



corrections were also given for unpolarized and polarized 
VCS observables both at low energies and in the Bjorken 
regime. The results will be an important tool for the analy- 
sis of present and forthcoming VCS experiments, in order 
to extract the nucleon structure information from these 
experiments. 



7 Conclusions and outlook 

It has been discussed how the real and virtual Compton 
scattering in different kinematical regimes provide new 
tools to extract nucleon structure information. 

It has been seen that for RCS at low energy, new accu- 
rate data have become available which not only allow the 
extraction of scalar polarizabilities of the proton, but also 
start to explore the spin polarizabilities of the nucleon. 
Those spin polarizabilities have been calculated recently 
to 0(p 4 ) in HBChPT. A fixed-t dispersion relation formal- 
ism was developed to extract the nucleon polarizabilities 
with a minimum of model dependence from both unpo- 
larized and polarized RCS data. The DR formalism was 
also used to obtain information on new higher order polar- 
izabilities of the proton, providing new nucleon structure 
observables and a new testing ground for the chiral calcu- 
lations. 

The VCS reaction at low photon energy maps out the 
spatial distribution of the polarization densities of the pro- 
ton, through generalized polarizabilities. Over the last few 
years, the VCS has become a mature field and a first ex- 
periment at MAMI at low energy has been successfully 
completed. In order to extract GP's from VCS data over 
a larger range of energies, a dispersion relation formalism 
is underway, providing a new tool to analyze such data. 
The DR formalism provides already results for 4 of the 6 
GP's, which can be confronted with chiral predictions. 

The RCS reaction at high energy and large momentum 
transfer is a tool to access information on the partonic 
structure of the nucleon. The PQCD predictions for wide 
angle real Compton scattering (90° in the cm.) show a 
strikingly different behavior than competing soft-overlap 
type mechanisms, and forthcoming experiments can teach 
us about the interplay of both mechanism at accessible 
energies. 

The VCS reaction in the Bjorken regime and associ- 
ated hard electroproduction reactions give access to new, 
generalized (skewed) parton distributions. The study of 



SPD's has opened up a whole new field in the study of 
nucleon structure. These observables unify two different 
fields as they interpolate between purely inclusive quanti- 
ties (parton distributions) on the one hand and between 
simple exclusive quantities (such as form factors) on the 
other hand. Besides the SPD's H and H, which reduce in 
the DIS limit to the quark distribution and quark helicity 
distribution respectively, there are two entirely new lead- 
ing twist SPD's (E and E), which cannot be accessed in 
DIS. The non-perturbative information contained in the 
SPD's is rather rich as they are functions of 3 different 
variables. In particular, the skewedness variable £ leads 
to different regions where one is sensitive either to quark 
distribution type information or to meson distribution am- 
plitude information in the nucleon. Through a sum rule, 
two of these SPD's (H and E) determine the quark orbital 
angular momentum contribution to the nucleon spin. An 
educated guess was shown for these SPD's, which was used 
to estimate the leading order DVCS amplitude. Further- 
more, the leading order meson electroproduction ampli- 
tudes were discussed and compared to the available data. 
In particular it was seen that in the intermediate W range 
(valence region) , a dominance of the handbag mechanism 
is predicted for p° electroproduction, which seems to be 
well confirmed by recent HERMES data. Furthermore, the 
extension of the formalism of the SPD's to the N — > A 
transition was discussed. The large N c limit allows to re- 
late the N ^ A SPD's to the N -> N SPD's. The trans- 
verse spin asymmetry was discussed as a promising ob- 
servable, likely to be less sensitive to higher twist effects. 

It is easy to foresee that the fields of real and virtual 
Compton scattering will show important activities in the 
near future both on the theoretical and experimental sides, 
and that an attempt to review them is very timely. It is 
hoped however, that the works initiated and discussed in 
this paper will stimulate further efforts on the theoretical 
and experimental sides to extend our knowledge of nucleon 
structure in new directions. 
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